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By W. Perl 


SUMMARY 

A simple approximate method is given for the calculation of 
isentropic irrotational flows past symmetrical airfoils, including 
mixed suhsonic -supersonic flows. The method is based on the choice 
of suitable values for the streamline curvature in the flow field 
and subsequent integration of the equations of motion. The method 
yields limiting solutions for potential flow. The effect of circu- 
lation is considered. 

A comparison of derived velocity distributions with existing 
results that are baaed on calculation to the third order in the 
thiclmess ratio indicates satisfactory agreement. The results are 
also presented in the form of a set of compressibility correction 
i*ules that lie between the Prandtl-Glauert rule and the von Karman- 
Tsien rule (approximately). The different rules correspond to 

different values of a local shape parameter '\/yCq, in which Y is 
the ordinate and C^ is the curvature at a point on an airfoil. 
Bodies of revolution, completely supersonic flows, and the signif- 
icance of the limiting solutions for potential flow are also briefly 
discussed. 


INTRODUCTION 

The problem of calculating compressible potential flows past 
aerodynamic shapes will be considered in this paper by the following 
method: An assumption is made as to the variation of the curvature 

of the streamlines in the flow field and the equations of motion, 
expressed in terms of the streamline curvature, are thereupon 
integrated . 

This basic method of calculating fluid flows is not new. It 
has been described in reference 1 for use in calculating pressure 
distributions on closely spaced airfoils in cascade. More recently 
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the method waa applied to the compreaaihle flow past an isolated 
airfoil and to the incompressible flow past a symmetrical airfoil 
in a closed channel. The results of reference 2 are compared with 
those of the present paper. 

The method is applied to isolated airfoils and the results 
are compared with those of references 3 and 4. The limiting solu- 
tion for potential flow by this method is identified with the 
"limiting line" phenomenon (reference 5) and discussed in relation 
to the flow through a converging -diverging channel. Application of 
the method to bodies of revolution is indicated. 


The flows calculated in this paper are of the steady, continuous, 
isentroplc, irrotational type. (See fig. 1.) The equations of 
motion are considered in the following form: 

Equation of irrotationality (reference 6, p. 43): 


Bernoulli's equation and equation of state for isentroplc flow: 


THEORY FOR SYMMETRICAL AIRFOILS 



( 1 ) 


Equation of continuity of mass flow: 


i 



( 2 ) 


1 


1 


P 



( 3 ) 


where 


V 


velocity at a point P of flow field 
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n distance measured along potential line from airfoil (streamline) 
to point P on stream].ine in f].ow field 

nQ perpendicular distance "between same two sti’ea]nline8 in free 
stream 

C curvature of streamline at point P, positive when streamline 

is convex in positive n direction. (C is the reciprocal 
of the radius of curvature.) 

p density 

p pressure 

7 ratio of specific heats 
M Mach number 


The si'.hscript o denotes free-stz'eam conditions. 

Consider the syimaetrical flow past a symmetrical airfoil section 
(fig. l). The potential lines of the flow pattern ai’e assumed 
straight and perpendicular to the free-stream direction^ or x-direction, 
thei’eby relaxing the condition of orthogonality between stream and 
potential lines. This assumption is exactly satisfied at the mid- 
chord station of the section if the section has fore-and-aft symmetry 
with respect to the mid-chord station. The assumption, in effect, 
renders the analysis for one choi-dwise station independent of that 
for another. 


Tile element of length dn in equations (l) and (2) is therefore 
replaced by the element of length dy in the y-direction and equa- 
tions (l) and (2) are written, respectively: 



dy = - 


ry 


dv 

Cv 

U) 

•£-^dy 
Po ^o 

(5) 


The differentials in equations (4) and (S) are understood to be 
taken in the y-dirocticn at constant x. The lower limits of inte- 
gration in equation (2) are on the streamline that coincides with the 
airfoil contour. The corresponding lower limits in equation (s) are 
therefore 0 and Y, respectively, whore I is the airfoil ordinate, 
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a function of the chordwlse location x. The upper limits of inte 
gration in equation (5) are on a streamline timt ■becomos parallel 
to the x-axis at infinity; that is: 

lin y = yo 

y— >oc 


Eenco^ the continuity condition (s)^ which for finite y^ Jq my 
he written 


fY ry 

n V 

fJo 

! dy = i 

.P 

dy - ! dy 

JO Jy 

Po "''’o 

JY 

becomes in the limit as y->cc 






1 

Y = I 

/ pv 

\ 

- i ; dy 

jy 

■.^1^0^ 0 

/ 


( 6 ) 


Vnthout loss of generality^ the free-stroam density Pq and 
velocity Vq are hereinafter censidored as uniter (or what is the 
sane, p and v are written in place of p/pQ and v/v^, respec- 
tively). Coii'hination of tho approxinato irrotational condition (4) 
and the approximate continuity condition (6) yields 



(pv-l) 

Cv 


dv 


( 7 ) 


in which the lower limit of integration, unity, is the free-stream 
velocity at y = oo, and tho upper limit V, corresponding to tho 
airfoil ordinate Y, is the unknora desired velocity at the airfoil. 


A streamline curvature fmiction is no^r to ho chosen. It must 
satisfy tho boundary conditions of known airfoil curvature Ca_ at 
tho surface of the airfoil and zero curvature at infinity » For 
convenience, in the integration of eqxiation (?), tho curvature C is 
chosen as a function not of the coordinates x, y of tho flow field 
directly hut of the velocity v(x, y). Tho function chosen is 


C = 


-'a 


Vv-i; 


( 8 ) 


in which tho airfoil curvature tho unkno\/n airfoil velocity 

V, and tho parameter t] (which will he discussed later) arc func- 
tions only of chordwise location x, hence are constant as far as 
tho integration in equation (7) is concerned. 
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The velocity v in the flow field is obtained by substitution 
of equaoion (8) into equation (4) and integration (at constant x)^ 


(V-1)^ n d v 


( 9 ) 


The streamline curvature function C(y) determined by equations (8) 
and (S) varies monotonically, as y increases at constant x, from 
the value at the airfoil where y = Y to the asymptotic value 

zero at y = oo. The manner of this variation depends on the known 
airfoil shape parameter YCg^ t2ae airfoil velocity Y, and the 
parameter t] . 


The airfoil velocity V is given in terns of the airf'oil shape 
parameter YCg^ the free-stream iiach number and the param- 

eter T) by equations (l), (3), and (s), which yield: 


YCa = (V-l)T 


^ V 1 - (v2-l)jY-l _ 

Jl v(v-l)’^ 


dv 


( 10 ) 


The parameter q is limited in its possible reinge of values by 
the conditions that must be satisfied infinitely far from the air- 
foil, that is, as v— ?1. These conditions are: 

(a) The curvat'ore C-*-0. 

(b) The distance y “—vOC • 

(c) The continuity integral (equation (lO)) is finite. Condi- 
tion (a) requires, by equation (8), that 0<q<cc. Condition (b) 
requires that the integral in equation (9) diverge as the upper 
limit V— >1, which it does for l<i'i<«x<. Condition (c) requires 
that the integral in equation (lO) converge at the lower limit, which 
it does for -a.<q<2. All three conditions therefore limit the 
permissible range of q to 


1< q< 2 (11) 

Equation (3.0) and condition (ll) represent the basic result of 
the present method. The velocity V at a point on the surface of a 
symmetrical airfoil at a given subsonic free-stroam liach number Mq 
is obtained by assuming that the Icnown data arc the shape param- 
eter YCa incompressible, or low-speod, velocity at 

that point on the airfoil. Tlio parameter q is first obtained from 
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equation (10) with iMq = 0 and 7 = It will appear from later 

applications that the value of T) thus determined falls within uhe 
range given hj expression (ll) over the region of the airfoil of 
greatest interest. With this value of p,, equation (lO) than gj.ves 
the velocity V at the point under consideration for the desired 
Mach number Mq. 

Tlie main uncertainty of the present method is represented hy 
the chosen curvature function given in equation (S). Tho approxi- 
mation involved in the use of this function has boon minimized "by 
tho method just descrihod of fixing the parameters V and q. By 
this method the streamline -curvature function satisfies an approxi- 
mate form of the equations of motion (equation (lO)) and yields the 
kno'vm exact value of the airfoil velocity in the incompressible case 
Furthermore, the curvature function exactly satisfies the boundary 
conditions of end values Cg^ and 0 and varies between these end 
val.ues in the correct general manner (for stations near the maximum 
velocity station), namely, monotonically. 


As a further condition on the chosen curvature function, the 
final solution given by equation (20) should reduce to the Prandtl- 
Glauert rule for siaali disturbance of tho free stream. Ihus, neg- 
lecting powers of v-1 and 7-1 equal to or higher than the second 
equation (lO) reduces to 


YCa = (7-1)"^ 


V !l - ^ (v+1)(v-1) - .... - 1 

L 6 ^ j 

v(v-l)^ 


dv 


T7 


= (Y-D'l 


(v-l)(l-Mo^) 

(v-1)^ 


QV 


r7 

= (l-Mio^)(7-l)’^ J (v-l)^"^ dv 


or, integrating and solving for 7-1, 

Equation (12) shows that, if q is adjusted to give the correct 
value of 7-1 for Mq = 0, or 


7i-l = ^/2-*q^/YCa 


(13) 
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then the compressible value for anj- subsonic Mq is given by 


v/l-Mo^ 


( 14 ) 


vhich is the Prandtl-Glauert rule. 

The foregoing considerations indicate that the results should 
not be critical with respect to choice of form of curvatui*e function. 
A somewhat different form of curvature function satisfying all the 
foregoing conditions is 


Comparative results based on this function will be discussed later. 

The basic relation (lO) connecting the velocity V (expressed 
as airfoil velocity increment V-l) at a point on a s;ymmetrical 
airfoil, the airfoil curvature parameter at the same point, 

and various free-stream Mach number Mq is shown graphically in 
figures 2, 3, and 4 I’or 7 = 1.4. The numerical data from which 
these curves were plotted are given in tables I and II. Included 
in tables I and II are some corresponding computations based on the 
curvature function (15). The integral in equation (lO) was evaluated 
by Simpson's rule and checked by the closed-form result obtainable 
in the case q = 1. The velocity increment above free-stream velocity 
V-l was plotted against hereinafter ca lle d the curvature 

parameter, rather than against YCa because */YCa, is proportional 
to the thickness ratio of the airfoil (see, for example, equation (19)) 
and V-l is therefore approximtely linear with respect to this 
quantity for sma.ll values (equation (12)). 

Figures 3(a) to 3(h) correspond to positive velocity increments 
above free-stream velocity, that is, evaluation of equation (lO) for 
V>1. Figures 4(a) to 4(h) correspond to negative velocity increments 
(V<1) such as produced on surfaces of negative cum^ature and ordinate 
(hence the negative sign attached to >JyCq^) . In this case the curva- 
ture fumetion, instead of the one given by equation (8), is properly 
taken as 


which allows equation (10) to bo evaluated without ambiguity as 
regards the berms to the power q . Tlie curves of negative velocity 



(15) 



8 


KA.CA ra Ko. 1328 


incroiaent (figs. 4(a) to 4(h)) may ho rogarcied as coritinuatlons of 
the corrospondlng curves of positive velocity increment (figs. 3(a) 
to 3(h) ) . 

It is of interest to note that a pai-amotor equivalent to the 

parameter ^A^Ca, was derived hy the authors of reference 8 from a 
dimensional consideration of equation (4) and used to correlate 
critical Mach number data for various symmetrical airfoil sections. 


APPLICATION AND EESUI.TS 
Kaplan Section 

Kaplan (reference 3) has calculated the compressible -flow zero- 
lift velocity distribution for a particular family of symmetrical 
aix’foils. The method used was an extension of the Ackeret method 
wherein the potential function is expressed as a power series in the 
thickness ratio of the section. The corresponding series for the 
velocity distribution was eva?.uated to the term in the third power 
of the thickness ratio. A limiting value of free -stream Mnch number 
was found^ for a given thicknoss ratio, above which the terms of the 
power series that were calculated (the first throe) indicated a 
probable failure of the series to converge. This free-stream Mach 
number was presumed to constitute an upper limit for the existence 
of a continuous potential flow. 

The method of the present analysis was tested by determining the 
velocity distributions for tho Kaplan section from figures 2, 5, 
and 4 (dra^m to a scale commensurate with the precision of the data 
of tables I and II). As an example of the procedure used, the 
velocity increment V-1 for the mid -chord location x = 0 of the 
Kaplan section of thickness ratio 0.10 was obtained as follows. The 
section ordinates and corresponding curvature parameter 
given in figure 5. Tho incompressible velocity distribution Vj_-1, 
obtained by conformal mapping, is shown in figure 6 (Mq = O). At 
X = 0 the values ^/YCa = 0.1925 and Vp-l = 0.1667 from figures 5 
and 6, respectively, correspond in figure 3(a) to an interpolated 
value of T) = 1.297. For these values of t) and \/YCa , the 
velocity increments V-1 for values of Mq of 0.5, 0.7, O.S, 0.85, 
and 0.9 were interpolated from figures 3(c), 3(e), 3(f), 5(g), 
and o(h), respectively. Velocity increments were obtained in this 
manner in tho chordwise range 0<x< 0.616. In tho chordwiso range 
0.616 <x< 1.0, the values of ^/YCa are indicated in figure 5 as 
imaginary, resulting from a positive ordinate and a negative curvature. 
The theory presented, based on the curvature function of equation (8), 
cannot handle such values of the curvature parameter . Approximate 
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Gompre3f3i'D.i.e velocity increments were obtained in this chordwise 
range by fii'st noting that the parameter t) increased toward the 
value 2.0 vath increase i.n chordwise distance from the center 
(fig. 5). Now the limiting value = 2.0 corresponds to the 
Prandtl-G-lauert rule, as will be shown voider COtiPRESSIBILITY COSEEC- 
TION RULES. In the range 0.616 <x< 1.0; therefore, the compress- 
ible velocity inci'eraents were calculated from the incompressible 
values by equation (14). 

The velocity distributions thus obtained for the Kaplan section 
of tliickness ratio 0.10 are shown in figure 6. The distributions for 
Mq = 0.85 and Mq = 0.9 do not extend all the way to the mid -chord 
location but come to an end (with infinite slope) at the chordwise 
locations x = 0.145 and 0.390, respectively. The immediate reason 
for this behavior is evident from figure 2. For example, at the 
limiting location x = 0.145 for Mq = 0.85, the value of q was 
1.312 (fig. 5). On the set of basic curves for = 1.312 and 
similar in appearance to figure 2, a vertical line drawn at 
v/yC^ = 0.181 corresponding to x = 0.145 would be tangent to the 
Mq = 0.85 curve, at W'hich point the velocity increment V-1 would 
be 0.550. No solution exists at this value of for Mq higher 

than 0.85; or, for fixed Mq = 0.35, no solution exists for higlier 
values of v^YC^ svich as correspond to chordwose locations closer to 
mid-chord than x = 0.145. The points of infinite slope on the basic 
curves of figure 2 coirrespond to a limiting solution for potential 
flow by the present method. This phenomenon, hereinafter called the 
"potential limit phenomenon," is discussed in appendix A. 

The potential limit points for the Kaplan 10-percent section at 
Mq = 0.85 and 0.9 were actually obtained from accurately determined 
plots of the infin5-te slope loci of figijre 3. These plots are given 
in figure 7. The vacuum-line boundary- curve in figure 5 corresponds 
to p/pQ =0 in equation (3). The intersection of the curve of q 
against \/YC^ for the Kaplan 10 -percent section (fig. 5), with, for 
example, the Mq = 0.85 contour in figure 7(a) determined the poten- 
tial limit values of q and -/YCg^ ; hence by figure 5 determined 
the chordwise location at which a potential limit point existed for 
I'Iq = 0.85. The potential limit value of q then determined the 
potential limit "v-1 by figviro 7(b). The free-stream ^fe.ch number 
at which the mid -chord location x = 0 is a potential limit point, 
that is, the lowest Mq at which a potential l.imit occurs, is 
indicated by point A in figure 7(a). By interpolation this limit 
value of Mq is estimated as 0.843 and the corresponding V-1 
(fig. 7(b)) as 0.573. For comparison, the lower critical Mach 
number (the lowest froe-stream Mach number at which sonic velocity 
occurs on the airfoil) was determined as 0.748. 
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The velocity distribution at the lowest potential limit Mach 
number, obtained by interpolation from the other velocity distribu- 
tions, is shown in figure 6. The velocity gradient at mid-chord 
appears to be finite and different from zero. The reason for this 
behavior can be seen from the expression for the velocity gradient 

dV BV BV du 

— = — ■==: + (16) 

dx bY yCg_ dx on dx 

in which the first parti al d erivative is taken at constant t| and 
the second at constant ^YCa. (The acceleration of the fluid along 

the surface a of the airfoil is V — . Hence, the following dis- 

ds 

cussion applies also to the fluid acceleration.) At a potential 
limit point both BV/o^YC^ and Bv/3n are ixifinite (fig. 2). If 

d^YC^/dx and dr)/dx are not zero at the ohordwise station corre- 
sponding to a potential limit x)oint, the velocity gradient dV/dx la 
infinite there, as is the case with the velocity distributions for 

Mq = C.85 and 0.9 in figure 6. If d-JlCg^/dx and du/cijc are zero, 
which is the case at the mid-chord station of the Kaplan 10-percent 
section, the occurrence of a potential limit at this point loads to 
an Indeterminate expression for the velocity gradient in equation (16). 
Closer analytical and graphical examination indicates the finite 
gradient shovm in figure 6 (hence, a finitely discontinuous change in 
fluid acceleration across the mid-chord station) . 


The lowest potential limit values of the corresponding 

values at mid-chord of V-1, and the local Mach number M for 
symmetrical sections of three thickness ratios (t = 0.05, 0.10, and 
0.20) were computed. The value of M was obtained by combining 
Bernoulli’s equation in the form of equation (3) and in the form 


resulting in 


y-1 9 

/ 1 + pi2 

^ 

^ (y2_l) 


(17) 


(18) 


The values are listed in the following table and compared with values 
obtained by Kaplan in reference 3: 
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Curvature method 

Kaplan luethod 

t 

Mo 

V-1 

M 

Mo 

V--1 

M 

0.05 

0.898 

0.324 

1.269 

0.890 

0.249 

1.164 

.10 

.843 

.573 

1.491 

.833 

.456 

1.320 

.20 

.763 

1.200 

2.257 

. 743 

.864 

1.625 


The limiting values of by the two methods agree fairly closely; 

the values of V-1^ hence also less so. The less satisfactory 

agreement of V-1 is to be expected from the rapidity with which 
V-1 varies with Mq in the neighborhood of the potential limit 
solution (fig. 6). 

The comparison of the velocity distributions for the 19-percent 
thiclmess section' by the two methods is given in figure 8. The 
velocity distributions for Mq = 0.75 (fig. 8(b)) and 0.83 (fig. 8(c)) 
were obtained by interpolation from cross plots of the velocity dis- 
tributions of figure 6, guj,ded by the potential limit points previ- 
ously determined. In the region of greatest interest on the contour, 
namely the supersonic region given approximately at Mq = 0.83 by 
0<x<0.4, the velocity by the Kaplan method increases with Mq at 
a greater rate than the velocity by the curvature method. A conven- 
ient criterion of the accuracy of velocity distributions in local 
supersonic regions has been pointed out by Tsien and Feeler, namely, 
if a velocity distribution indicates a local supersonic region, a 
velocity distribution can be derived in this region that must be 
greater than the original velocity dlstribiition; the difference 
between the two distributions decreasing as the extent of the local 
supersonic region increases. This greater velocity distribution 
for the supersonic region is the well-known Prandtl-Meyer solution 
for the flow over a curved surface (reference 8(a)). The Prandtl- 
Meyer velocity distribution extends, in the case considered here, 
from the chordwise location for a local Mach number M = 1 to the 
mid-chord location. The Prandtl-Meyer solution is obtained from 
the change in slope of the airfoil surface A® from the M = 1 
location to the point under consideration. This change in slope 
expressed in degrees is equivalent to the pressure number P, which 
is a ilmction of the local (supersonic) Mach number M. From a plot 
of slope 9 of airfoil surface against x, included in figure 5, 
and a table of P against Mach angle, it is therefore a simple 
matter to obtain the Prandtl-Moyer local Mach number distribution in 
the supersonic region. The local Mach number M and the free -stream 
Mach number Mq then determine the local velocity V by equa- 
tion (18) . If the pressure number P is defined in terms of the 
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flow deflection angle A0 as P = 1090 - the pressure numher 
inay be computed by the methods of reference 8 (a). The variation of 
pressure number and flow deflection angle with Mach angle is shown 
in the following table : 


PEESSUR2 rroi-IBEE AND MACH ANGLE IN PEANDTL -MEYER 
SOLUTION FOE SUPERSONIC FLOW 
[7 = 1 . 4 ] 


Pressure 
number, P 

Flow deflec- 
tion angle, 

A0 = 1000 - P 
(d.eg) 

Mach angle 

u = sin"^ i 
M 

1000 

0 

90° 00' 

999 

1 

670 43' 

998 

2 

620 00' 

997 

3 

580 XO' 

996 

4 

440 12' 

995 

5 

520 43' 

994 

6 

500 36' 

993 

7 

480 43* 

992 

8 

470 03' 

991 

9 

450 32' 

990 

10 

440 09' 

989 

11 

420 51' 

988 

12 

410 39 . 5 ' 

987 

13 

400 32.5' 


The Prandtl-Meyer solutions obtained in this manner for the 
velocity distributions by the Kaplan method and by the curvature 
method are shown in figure 8 (c). The velocity distribution obtained 
by the curvature method is evidently closer to the true distribution 
than that of the Kaplan method in the local supersonic region for 
Mq = 0.83, because the curvature velocity distribution is less and 
the Kaplan velocity distribution is greater than the corresponding 
Pz’andtl -Meyer distribution. Application of the Prandtl-Meyer solution 
to the potential limit velocity distribution by the curvature method 
for Mo = 0.843 (fig. 8 (d)) indicates a rapid decrease in validity 
of the results by the curvature method for free -stream Mach numbers 
close to the lowest potential limit Mach number. 
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From these and similar calculations for thickness ratio of 
5 and 20 percent, it is concluded that the curvature method gives 
results of at least the same order of accuracy as the Kaplan method 
over approximately the middle half of the Kaplan section. Over the 
rest of the section the Kaplan method may well he more accurate . 


The method of curvature may he expected to Increase in vali dity 
with decrease in the variation oi"* the curvature parameter 
along the airfoil. The symmetrical biconvex section (one formed of 
two circular arcs) was considered as an example which is more favor- 
able in this respect than the Kaplan section. The biconvex section 
of 0.20 thickness ratio was analyzed for its incompi’essible ve loc ity 
distribution by conformal mapping* The curvature parameter 
was calculated fTOm the relation 


in which Ymax is the section ordinate at the mid-chord location 
and t is the thickness ratio. The sect ion ordinates Y, the 
slope 0, the curvature parameter '\jYCa, and the parameter t), 
detennined as for the Kaplan 0.10 thickness ratio section, are 
given in figure 9. The velocity distributions for Mq = 0, 0.5, 

0.7, and 0.0 ai’e given in figure 10. The cui’vature method cannot 
handle velocities less than the free-stream value (unity), which 
correspond to positive values of the curvature parameter; these 
velocities were obtained by the Prandtl-Glauert rule. 

The lowest free-stream Mach number Mq for which a potential 
limit occurs at mid -chord is 0.790. The corresponding values of 
local velocity increment and local Mach number are 1.819 and 1.704, 
respectively. The Prandtl-Meyer solution (fig. 10) indicates 
inadequacy of the curvature method in this case at a free-stream 
Mach number somewhat less than 0.8. The local velocity increments 
of the biconvex section of thickness ratio 0.20 are higher than 
those of the Kaplan section of thicloiess ratio 0.10. The ratio of 
the increments of the biconvex section to those of the Kaplan section 
of thickness ratio 0.10 are, in general, less than the ratio of the 
thicknesses of the two sections. 


Biconvex Section 



CONDITIONS IN FIELD OF AIRFOIL 


The variation of streamline curvature and of local velocity 
with distance from the airfoil is given by equations (8) and (9). 
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For the purpose of illustration, it will be sufficient first to 
make the simpler calculations corresponding to the curvature func- 
tion (15) in place of (6). Sguations (4) and (15) yield 


which, after integration and adjustment of the constant of integra- 
tion to satisfy the boundary conditions, becomes 


The curvature variation is given by eq,uations (15) and (20) . For 
= 1 these equations are indetejcminate . In this case, equa- 
tions (4) and (15) (or the well-known limiting form of equation (20)) 
yield 


The curvature and velocity variations given by equations (15), 
(20), and (21) are shown in figure 11 and 12, respectively, for 
various values of q and Mq. The value of YCg, of 0.03704 chosen 
for the calculation was that corresponding to the mid-chord location 
of the 10-percent thick Kaplan section. The corresponding values of 
V for the various were those previously computed for the Kaplan 

10-percent thick section. Although these values of V corresijond 
to only a single value of t| based on the curvature function (8), 
the use of these same values for various r| and with the curvature 
function (15) give the trend of curvature and velocity variations 
sufficiently well for the purpose of illustration. The curvature 
and velocity variations of figures 11 and 12 show the expected 
trend with free -stream Mach muaber Mq, namely, a slower decrease 
to free -stream conditions as Mq is increased. The limitation of 
a single curvature function such as equation (6) or (15) is apparent 
from figure 13, for it yields basically the same kind of curvature 
variation in the supersonic region as in the subsonic region. A 
more rigorous analysis should take into explicit account the different 
typo of curvature variation found in supersonic flows. 

The extent of the local supersonic regloias in the field of the 
10-pei’cent thick Kaplan section was next calculated by the more 
appro'priate equation (8), using the data derived in connection with 


dy = - 


Uoge V)^ d(loge v) 

.logg v)’^ 




( 20 ) 



Ti = 1 (21) 
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the velocity distrlhutions. The results are shown in figure 14^ in 
which the local Mach number M = 1 boundary is plotted for Mq = 0.75, 
0.8, 0,83, and 0.843 in terms of airfoil aemichord as unit distance. 

The rapid increase in lateral extent of the supersonic region 
with increase of free -stream Mach number, particularly near the 
potential limit, is apparent in figure 14. At the potential limit 
Mq = 0.843, the local supersonic region extends laterally into the 
flow field a distance of about 1 chord. 

The finite nonzero slope at x = 0 of the sonic boundary at the 
potential limit should also bo noted. This finite slope is associated 
with the finite velocity gradient at the 3\i.rface of the Kaplan section 
previously discussed. It appears from equation (9) that the slope of 
a constant velocity (hence, for given Mq, constant M) boundary 
will be associated with the surface velocity gradient at the same 
value of X in being finite or infinite. Hence, if the lowest poten- 
tial limi t sol ution occurs at a i>oint on a symmetrical airfoil at 
which d /\/yCa,/dx, dq/dx are not zero, then the consequent infinite 
velocity gradient at that point on the surface will cause a cusp in 
the M = 1 boundary in the field. This cusp (the possibility of which 
was suggested by L. Richard Turner of the NACA Cleveland staff) amounts 
to an envelope of the Mach lines in the (supersonic) neighborhood. 

This behavior is in agreement with a known property of potential limit 
solutions, namely, that a potential limit point in a flow field lies 
on an envelope of Mach lines (reference 5) . 


CIRCULATOEI FLOW PAST CIRCULAR ARC MEAN CAMBER LINE 

The curvatvare method was next applied to the calculation of the 
type of velocity distribution that produces "design" lift, that is, 
lift without a velocity peak near the nose of the airfoil. The 
circular arc camber line at zero angle of attack was chosen for this 
calculation in order that a comparison analogous to that for the 
Kaplan section could be made (reference 4) . The camber ratio assumed 
was 4 percent, corresponding to an incompressible design lift coef- 
ficient of 0.520. 

A difficArlty of principle arises in obtaining the curvature 
parameter VYCa for the circular arc section. The ordinate Y as 
developed in the curvature method is actually the component in the 
y-dlrection of the distance between a point on the airfoil contour 
and a point on the airfoil streamline infinitely far from the airfoil 
(the airfoil streamline is the streamline that includes the airfoil 
contour). In the case of the symmetrical airfoil, this projected 
distance is the airfoil ordinate as measured from the chord line. 

When lift is produced, however, a point on the airfoil streamline 
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infinitely far from the airfoil in the x-direction is also infinitely 
far from the airfoil in the y-direction. Nevertheless, Inasmuch as 
the slope dy/dx of the airfoil streamline for the circular arc 
rapidly approaches zero with increasing distance from the circular 
arc, it might be expected tlmt the velocity distribution on the 
circular arc would not be critical with respect to the finitely 
distant point on the airfoil streamline from which the ordinate Y 
is measured. 

This situation was studied by obtaining in the usual manner the 
compressible velo city increment at the mid-chord location, for curva- 
ture parameters y^YC^ corresponding to various values of Y as 
determined by points on the incompressible (Mq = 0) airfoil stream- 
line at various distances from the circular arc. The resulting 
velocity increments (V-l)iuj^ are plotted for various free -stream 
Mach numbers Mq in figure 15 against a quantity f characterizing 
the distance along the airfoil streamline . The quantity f is the 
ratio in incompressible flow of the velocity decrement at a point 
on the airfoil streamline to the maximum velocity increment on the 
circular arc. As the variable point on the airfoil streamline 
approaches the circular arc, the value of u decreases from 1.28 at 
f = 0.0032 to the lower limiting value 1 at f = 0.0102 correspond- 
ing to a distance of 0.67 chord from the extremity of the circular 
arc. At points on the streamline closer to the circular arc there 
is no valixe of n that yields the knoTO incompressib^ maximum 
velocity increment for the corresponding value of VYC^. Compressible 
maximum velocity increments for these points were thereupon arbitrarily 
obtained from figure 2 using p = 1. 

Included in figure 15 are the maximum-velocity increments 
calculated from the formulas and constants given in reference 4. 

The maximum-velocity increment by the curvature method is seen to 
be always less than Kaplan's value and not to vary greatly with f. 

Guided by the results for the Kaplan 10-percent thick symmetrical 
section, which had about the same incompressible maximum velocity 
increment as the 4 -percent camber circular arc mean line, and also 
by the comp aris on with Kaplan's results in figure 15, the cxirvature 
parameter -/y^ obtaining complete velocity distributions was 

determined with respect to a reference point on the airfoil streamline 
at which the incompressible velocity decrement was 1 percent of the 
maximum velocity increment (f = 0.01, 0.68 chord from leading edge). 
The basic data are given in figure 16 and the resulting velocity 
distributions are compared with the corresponding results by the 
Kaplan method in figure 17. Tlie velocity distributions at Mq = 0.8 
by both methods are leas than the Prandtl-Meyer solution by about 
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the samo amount and no definite conclusion as to the comparative 
accuracy of the tvo methods is indicated. The uncertainty introduced 
into the res ult s hy the amhiguity of determination of the cxi'vature 
parameter V^Ca could presumably he further resolved hy additional 
comparisons with results hy other methods or hy calculation of higher 
approximations hy use of the streamline curvature. 


A compressihi.lity -correction rule may he defined as a rule hy 
which the velocity or pressure at a point of a compreasihle flow field 
is expressed as a function only of the low-speed, or incompressihle, 
velocity or pressure at that point and the free-stream Mach number. 
Thus, the curves of figure 2 when cross-plotted against Mq constitute 
a set of compresaihility-corroction rules for the velocity V. In 
general, each value of p yields a compressibility -correction rule. 

The rules for p = 1.0 and p = 1.8 are shown in figure 18 in terms 
of the pressure coefficient C^, defined as 

XT 


in which v/Vq is given in terms of the velocity V by equation (3). 
Also shown in figure 18 are the potential limit curves, obtained from 
figure 7(d), the sonic line (M =1 in equation (17)), and the abso- 
lute limit or vacuum line (p/P© = ^ equation (22)). The differ- 
ence between the rules for different values of ^p illustrate the 

allowance for differences in airfoil shape (y’YCa.), which correspond 
to the same low-speed pressure coefficient. 

Each compressibility-correction rule has its own potential limit 
curve. The limit curve for p = 1.3 intersects the vacuum line and 
would extend to local velocities higher than that corresponding to 
zero pressure, which of course is impossible. Hence, where a potential 
limit curve intersects the vacuum line, the vacuum line becomes the 
limit curve. 

In figure 19 comparison of the compressibility correction rules 
by the curvature method is made with the rules of 

(a) Prandtl-Glaucrt, equation (14) 

(b) Kaplan (reference 3) results for mid-chord location 


COMPRESSIBILITY -COKRECTION RULES 


- 1 

P-Po Po 


( 22 ) 



(c) The rule for p = 1.0 based on function (15), C = Cg. 
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(d) Greene (reference 2) 

(e) Garrick-Kaplan^ "arithmetic -mean" laile of reference 10 

(f) von Karman-Tsien_, equation (62) of reference 9 

Portions of the potential limit curve corresponding to each rule 
are indicated . 

The compressibility-correction curve corresponding to the 
curvature function (15) rises somewhat more steeply than that corre- 
sponding to equation (8) for the same value of r, . It should be 
noted that the two compressibility correction rules for n = 1 do 
not indicate the difference in cal culat ed pressure coei-ficient at a 
given point on an airfoil (given '\jYCg_) as a result of choice of 
curvature function. In figure 19 the curve for curvature func- 
tion (15) corresponds to = 0.305 and that f or c urvature func- 
tion (8) to VYC^ = 0.317. For a given value of and low-speed 

pressure coefficient, the results by both curvature functions differ 
negligibly. 

The Prandtl-Glauert rule has been drawn in figure 19 to the 
vacuum line (M = oo) because this rule can be regarded as the limit 
of the rules derivable by the present method as r, approaches the 
value 2.0, This fact becomes evident from equation (10), in which, 
as T] approaches 2.0, the velocity V must approach unity in order 
for the integral to converge. The approximation leading to the 
Prandtl-Glauert rule (eqi;ation (14)) can therefore be made. As q 
approaches 2.0 the vacuum line becomes a greater and greater portion 
of the complete potential limit curve (fig. 18) until in the limit 
the vacuum line becomes the entire potential limit curve. (See 
figs. 7(b) and 7(d).) 

The Greene rule, on the other hand, can be regarded as the 
limiting rule obtained by alloi/ing q to approach zero. For q = 0, 
equation (10) becomes 


This last equality, namely V as a function of and Mq, is 

Greene's rule. As has been shown, the significance of q = 0 is 
that the curvature C becomes zero at a finite distance from the 
airfoil and remains zero for greater distances. Although this 
. circumstance might imply a severe limitation on the validity of 


1 



1^1 - 2:^ Mq2 (v2-1)1 


dv - loge V (23) 


= (Vi - 1) - logg Vi 
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Greene's rule, this particular derivation is actually more restric- 
tive than necessary for 3t leads not only to Greene's inile but to 
the additional equation involving YCg,. A more general discussion 
of Greene's rule is given in appendix B. 

In general, the compressibility-correction curves derived in 
this paper lie between the Prandtl-Glauert curves on the one Imnd 
and (apriroximately) the von Karman-Tsien curves on the other. It 
may be emphasized here that the presentation of the present results 
in the form of a set of compressibility correction niles does not 
im]ply an equivalence of these results to a simple speed distortion 
of the flow field in going from incompressible to compressible flows. 
Evidently any derived set of compressible -flow patterns for various 
free -stream Mach numbers can be compared with the corresponding 
incompressible -flow pattern by means of a set of compressibility 
correction rules. A simple speed distortion implies the existence 
of only a single compressibility correction rule from which 
compressible -flow boundary velocities are obtained from given 
incompressible -flow boundary velocities regardless of the shape of 
boundary tnat prodi’ces the incompressib.le-f3.ov velocities. The 
present results, however, yield different compresoibie-flow veloci- 
ties for the sam_e incompressible velocities depending on the shape 
of airfoil (^YCa) that produces the incompressible velocities. 


SUPERSOIIIC FLOW 

In completely supersonic isentropic potential flow the same 
equations of motion hold as in the subsonic case, namely, in the 
form assumed in this paper, equations (l), (2), and (3). The appli- 
cation of these equations, however, to the calculation of the super- 
sonic velocity distribution on airfoils Is in cor';ain respects dif- 
ferent f>.’om the subsonic calculation. The differences as well as 
similarities will be illustrated by derivaiion of the linearized 
(small perturbation) equation for the velocity distribution on an 
arbitrary thin airfoil in supersonic potential flow. 

For small perturbation of the free stream the approximation for 
the density ratio that led to equation (12) is 

pv - 1 = (1 -Mq^)(v-1) 

Equation (7) therefore becomes 

PY 

Y = -(Mo^-l) / dv 


(24) 

f 


u 


1 


(25) 
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The curvature C is by definition dp/ds where s is the distance 
along a streamline and P is the (positive clockwise) angle between 
the streamline and a fixed (say the free-stream) direction. Equa- 
tion (25) may therefore be written 


2-1) / (v-1) ^ ds (26) 

Jo 

in which the path of integration is determined as follows: The 

assumption is made that at a point in the flow field any field 
quantity, such as v, dv/dp, or ds, is constant along a straight 
line connecting the point with a point on the airfoil and making an 
angle n with the free-stream direction. (See sketch.) 




Completely supersonic potential flow 

By this assumption a correspondence is set up between points in the 
field and points on the airfoil for which the integrand in equa- 
tion (26) has the same value. Hence the original path of integra- 
tion from infinity to a point on the airfoil in a direction normal 
to the airfoil can be replaced by a path along the airfoil itself, 
fix>m the leading edge to the same point; at a distance s^ from the 
leading edge (the contribution to the integral from points upstream 
of the leading edge is zero because for such points the straight 
line at angle p on which the field quantities are constant does 
not intersect the airfoil and extends infinitely far ahead of the 
airfoil, where free-stream conditi.ons exist). 
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Differentiating "both sides of equation (26) with respect to the 
distance Sa along the airfoil contour yields 


dY 

isa 




But -I- 
ds 


= - P- 


Hence 


a 


v-1 dv _ 1 

^ Mc.^-1 


(27) 


(28) 


Inasmuch as all variable s have been asanmed to differ from their 

V ~1 

free -stream values by small amounts, the quotient — ^ can be set 
equal to the derivative dv/dp. Equation (28) therefore becomes 


^ = ± 1 


(29) 


Finally the significance of the angle p is obtained from equa- 
tions (4) and (29), and the sketch of completely supersonic potential 
flow, 


tan p = ^ = - ^ 

ds Cvds dp 

. ± -pi= (30) 

or 

sin u = (31) 

Hence n is the Mach angle and equation (29) can be written 

~ = ± tan |i (32) 

dp 

which (recalling that the free -stream velocity is unity) is the 
basic solution for linearized two-dimensional supersonic flow 
(reference 8(b)). All the properties of linearized two-dimensional 
supersonic flow can be deduced from the preceding derivation (these 
properties were, of course, a guide in setting up the derivation) . 
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The fact that the lineari 2 ed theory for completely supersonic 
f?.ows relates streamline curvatures at points connected by Mach 
lineo (characteristics) would Indicate that in obtaining higher 
approzi mat ions to the flow pattern in large local or completely 
supersonic regions the curvature function should properly be spec- 
ified along Mach lines rather than along normals to the free -stream 
direction. This circumstance indicates a possible refinement of the 
present treatment for mixed subsonic -supersonic flows. Consider a 
mixed subsonic -supersonic flow pattern calculated by the curvature 
method. The curvature in the local supersonic region can be assumed 
constant (and equal to the corresponding airfoil curvature) along 
the characteristics emariating from the airfoil at the appropriate 
local Mach angle. The curvature function thus determined is contin- 
uously joined to the curvature function previously calculated at 
the boundary of the supersonic region. As a guide in this process 
the curvature of the streaimlines determined by the original calculation 
could be used (this information would also be of use in determining 
higher approximations tliroughout the entire flow field) , The equa- 
tions of motion (1), (2), and (3) are thereupon integrated, graph- 
ically or otherwise. If the local supersonic region is small, this 
procedure may not yield more accurate results than the original 
calculation. It may, however, yield a closer evaluation of the 
accuracy of the original calculation than the Prandtl -Meyer solution 
applied to the local supersonic region. 

The preceding remarks also make evident the possible occurrence ' 
of a curvatvire maximum away from the airfoil (in a direction normal 
to the free stream) . Such a maximum will probably occur in a large 
local supersonic region when the point under consideration on the 
airfoil is preceded along the airfoil by points of greater airfoil 
curvature . 


The method given for the two-dimensional flow can be applied 
similarly to the case of axially symmotric flow over a body of 
revolution. Thus, if the s,'\Tnmetrical section (fig. 1) is considered 
as a meridian section of a body of revolution the contiiuiity con- 
dition (equation (6)) becomes (free-stream velocity and density are 
unity). 


BODIES OF KEVCLUTION 



pco 


(pv-1) y dy 


(33) 
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The irrotationallty condition (equation (1)) and the Bernoulli -state 
equation (3) reine,in unchanged. If the curvature function (8) is 
assumed, substitution of equations (3), (4), and (9) into equa- 
tion (33) yields an equation analogous to (10). From this equation 
a set of basic curves analogous to those of figure 2 could be derived 
by computation. (The computations would be lessened by use of the 
curvature function (15) instead of (8) because equation (20) is 
easier to evaluate than equation (9).) 

It ■will be sufficient here to Indicate the results in the case 
of small disturbance of the free stream, or v — #1. In this case 
equations (9) or (20) yield 


Substitution of equations (24), (34), (4), and (8) into (33) gives 


It may be shown as before that the parameter q is now restricted 
to the range 


The parameter tj can again be taken as the value that in equa- 
tion (35) yields the exact (known) velocity increment for Mq = 0. 

Equation (35) differs in form from the corresponding expression 
for two-dimensional flow given in equation (12) by virtue of the 
second teim in the brackets. The effect of this term is to reduce 
the velocity increment for a given Mq, that is, to reduce the 
effect of STjeed on the local velocities and pressures. For example, 
the increase with free -stream Mach number Mq of the maximum velocity 
increment on slender ellipsoids of revolution is by equation (35) 
about 70 percent of the increase given by the Prandtl-Glauert rule. 

If the curvature function is chosen as 



(34) 



(35) 


1 ^ Ti < 1.5 


(36) 


J (1)^ 


(37) 
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then with | = 1 a Dacic equation the same as (10) is obtained^ 

yc 

with the exception that — 2: replaces YC^ on. the left-hand side. 


Hence the same ba-sic computations as already made for the two- 
dimensional case could be used here (with 1 = p < 2) . Although a 
more gradual rise of velocity increment with Mach number is thus 
obtained than for the same section in two-d:Unensional flow, the rise 
is greater than given by the Prandtl-Glauert rule. The Prandtl- 
Glauert rule is, in fact, again the limiting rule for small disturb- 
ances. It seems, however, that the Prandtl-G.lauert rule may over- 
estimate the effect of subsonic compressibility speeds on slender 
bodies of revolution. 


Perhaps the most reliable way of obtaining compressible velocity 
distributions for bodies of revolution by use of a single curvature 
function is to use the function (37) with q and | adjusted in 
each case to satisfy both the knoAvn incompressible value and the 
compressible value for infinitesimally smq.ll disturbance of the 
stream, the compressible va.lue being considered as known or obtainable 
from the general liiiear-perturbatlon theory of compressible fluids. 

The permissible values of t] and i that satisfy the required 
conditions at infinity lie in the acute-angled sectors of the | 
against t| plane bounded by the lines p * 1 and | = 2t) - 3. 


SUlvDiMlY OF EESULTS 

The present study of compressible potential flows past aero- 
dynamic shapes indicates the following: 

1. The method presented for the calculation of compressible -flow 
velocity distributions yields results for sjnimetrical sections in 
satisfactory agi’eement with existing results baaed on calculation to 
the third order in the thickness ratio. 

2 . The resu.lts can be presented in the form of a set of 
compressibility-correction rules that lie between the Prandtl- 
Glauert rule and the von Karman-Ts.ien rule (approximately). The 
different rules correspond to different values of a local shape 
parameter ^\/YGa, in which Y is the ordinate and Ca is the 
curvature at a point on an airfoil. 

3. The effect of circulation at design lift conditions, that 
is, without velocity peaks, can be taken into account. 

4. Conditions in the field of the airfoil can be calculated 
simply. 
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5, The method given for two-dimensional flow can be applied 
also to bodies of revolution. 

6. The general method of using the streamline curvature appears 
applicable to any subsonic or supersonic flow problem in which a 
satisfactorily accurate estimate of the curvature of the streamlines 
can be made in the portion of the flow field of interest. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, May 24, 1946. 
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APPSI'iDIX A 


TEE POTENTIAL T.WP PHENOMENON AND THE ANALOGY 
WITH FLOV/ THROUGH A CHANNEL 

The points of infinite slope on the hasic cusrves (figs. 2 
and 3) have heen seen to correspond to limiting solutions for 
potential flow by the present curvature method. These points cor- 
respond to points on an airfoil at which the local Mach number M 
is greater than 1 and at which the fluid acceleration is infinite 
if dv/YCg^/ds and dii/ds are not zero. Furthermore the M = 1 
boundary in the flow field contains a cusp if the fluid at the 
laterally corresponding point on the airfoil has infinite accel- 
eration. This cusp constitutes an envelope of the Mach lines in 
the (supersonic) neighborhood. These properties permit identifi- 
cation of the points in the flow field corresponding to points of 
infinite slope on the basic curves >/ith points on the limiting line 
of reference 5. 

In the detearmination of the essential reason for the existence 
of a limiting potential flow solution, as well as the significance 
of the possibilities (VYCa/'is = di^/ds = 0, the analogy with the 
flow through a converging diverging channel is illuminating. Con- 
sider first the one -dimensional flow through the channel. The 
equation of continuity can be ■t/ritton 

pVA = 1 (38) 

in which A is the cross-sectional area of the channel and all 
quantities are expressed as fractions of their values at a reference 
station 0 upstream of the minimum section, called tho channel frec- 
stream station. Equation (38) yields with equation (3) a family 
of curves of velocity T against area A with channel free-stream 
Mach number Mq as parameter (fig. 20). Tiieso curves exhibit 
noints of infinite slope analogous to those of figure 2. If for 
the moment only subsonic flow in the channel is considered (branch 
BC in fig. 20), there exists at C a minimum channel area for , 
given Mq or for given channel area a raaximum Mq for continuous 
one -dimensional flow, Tho fluid acceleration along tho axis x of 
the channel is 
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in vliicli, ty squation (38), (l?)^ and (3) 

^ ^ Y__ 

A (m2-1) 


(40) 


At a point of infinite slope in figure 20, dV/dA = »; hence the 
fluid acceleration is infinite unless possibly dA/dx = 0. The 
limiting solu.tion occurs, by eqr.ation (40), at a local Mach number 
of unity. Thus, by a one -dimensional argument independent of the 
irrotationality condition a limiting solution analogous to that for 
the isolated airfoil has been derived. It may be noted that the 
upper -branch solutions CD (fig. 20), as well as the upper -branch 
solutions of the basic curves (fig. 2), correspond to the over- 
lapping supex’sonic flow patterns that have been obtained by the 
hodograph method. (See reference 11.) 

Hie analogy between the channel and the Isolated airfoil can 
be made still closer by consideration of the two-dimensional features 
of the flow in the minimum section of the channel. As the channel 
free-stream Mach number !>do is increased by increasing the over-all 
pressure difference across the channel through external means, the 
maximum local Mach number at the wall increases at a greater rate. 

I’or a sufficiently high but subsonic Mq, local supersonic regions 
appear in the neighborhood of the walls (fig. 21). This flow 
pattern correspoiids to the solution symmetrical with respect to the 
y-axis studied by Meyer (reference 12), G. I. Taylor (reference 13), 
and others. It is analogovis to the continuous mixed subsonic- 
sutpersonic flow pattern for isolated airfoils. The one -dimensional 
continuity treatment indicated an upper subsonic limit for Mq. 

The two-dimensional flow pattern for channels indicates a similar 
uppei’ limit on Mq and, in addition, provides the desired insight 
into the isolated airfoil case. 

The limit on comes about because the flow at station 0 

(fig. 21) must pass through the minimum section. The local mass 
flow Intensity pv is, however, a maximum at a local Mach number 
of unity and will be less than the flow intensitj'" at the channel 
free -stream station in portions of the local supersonic regions A3 
and FG, which increase as Mq is increased. When a further increase 
in extent of the local supersonic regions would result in a decreased 
mass flow through the minimum section from the cause Just indicated, 
then the mass flow has reached its maximum possible value and the 
channel is said to bo "chohed." 

The explanation for the limiting solution in the isolated air- 
foil case is similarly formulated. Along the potential lino AD 
traversing the local supersonic rogion (fig. l) the mass flow 
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intensity ratio — is a maximum at the sonic point C and less 
Po^o 

than unity over a section AB in the supersonic region. When the 
mass flow intensity, integrated over BD, produces an Insufficient 
mass flow to counterbalance the decreased contribution to the mass 
flow across AB in accordance with the req.uirements of continuity 
and isentropy (equation (6)), then the limit solution for point A 
on the airfoil has been reached. The flow field, though infinite 
in extent, can under this condition be said to be choked. Evidently 
a local supersonic region must exist before the limiting velocity is 
reached . 

The essentially one -dimensional continuity argument Just given 
is not entirely sufficient to prove the existence of a limit solution 
in the two-dimensional case. Equation (6) alone, for example, could 
always be satisfied by a suitable choice of v(y); or large stream- 
line curvatures resulting in large flow deflections might set in at 
supercritical speeds, thus destroying the validity of equation (6). 

These possibilities are eliminated by the irrotationality condition, 
which controls and limits the lateral variation of velocity in the 
flow field. Thus, the isolated airfoil limiting solution is actually 
pi’oduced bj" the combination of irrotationality, continuity, Bernoulli’s 
equation, equation of state, and boundary conditions in equation (10). 

The condition of completely irrotational flow is not, however, abso- 
lutely necessary for a limit solution. Eotational flows characterized, 
for example, by almost any function of v on the right-hand side of 
equation (1) would yield limit solutions, at least for small values 
of the function (rotation) . The equations of motion in the form 
analyzed appear, in fact, to offer a convenient means of including 
rotational effects and effects of changes in the equation of state. 

Finally, some remarks are made concerning the possible relation of 
the actual shock wave on an airfoil in the mixed subsonic -supersonic 
(supercritical) flow regime to the potential limit solution. As noted 
in references 5 and 14, the observed shock on an airfoil in the super- 
critical flew regime appears to be formed as a result of the ever-present 
random pressure disturbances, some of which travel upstream at relative 
sonic speed and pile up to form the downstream boundary, roughly speaking, 
of the local supersonic region. In other woi-ds, shock would not arise 
in a completely steady flow; and an increase of free -stream Mach number 
(a nonsteady effect) would presumably permit attainment of the potential- 
limit solution. Although it seems possible that the potential-limit 
solution might in some cases limit the local Mach number at which the 
shock stabilises, the maximum local Mach numbers corresponding to the 
potential limits of this paper are for the most pai't greater on normal 
airfoils than those at which the shock stabilizes. Hence, the actual 
shock normally prevents the potential -limit solution from being reached. 
There appears to be no direct relation between the two phenomena. 
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Th3 question as to what would happen if the potential limit 
solution could he attained in some manner on an isolated airfoil 
can peifiaps best he ans’.rerod hy again considering the convergjng- 
diverging channel. In the channel the limit solution can not only 
he reached (point C in fig. 20) hut the second solution (branch CD) 
can also exist. The mechanism that produces this second solution is 
the over-all pressure difference across the channel^ applied through 
external means. This p?.-es 3 ure difference, when sufficiently increased 
causes the shock, vrhich has fomied in the local supersonic regions 
in a manner similar to that on an isolated airfoil, to move dcTO- 
stream from the minimum section as a more or leas nonml shock 
spanning the channel. The region of the channel between the minimum 
section and the shock contains the second sol'ution. Thus, the 
upper branch solution CD can exist as a continuation of the lower 
branch solution BC with no shock in the neighborhood of the limit 
solution, point G. The fluid acceleration of eqi’.ation (39) is 
finite at the limit point which occurs at the minimum section 
dA/dx =0. In the isolated airfoil case the analogous occurrence of 


a potential limit solution at a point where 


<V^^a _ dri _ 
ds " ds ■ 


(equation (25)) leads to the conjecture whothei’ a similar transition 
through the potential limit solution could not be effected at such a 
point. In normal isolated airfoil operation the only mechanism 
available for increasing the local Mach number is to move the airfoil 
faster. The over -all ambient pressure remains atmosphei'ic. The fact 
that in the channel the necessarily asymmetrical boundary conditions 
of over -all pressure difference (t]ie same channel area at beginning 
and end of the channel is assumed) can produce an asymmetrical flow 
pattern, containing a transition tlirough the limit solution would 
therefore indicate that on an isolated airfoil subject to the 
symmetrical boundary condition of constant atmospheric pressure a 
transition through the potential limit solution could not be effected. 
Artificial moans, however, such as the proper combination of airfoil 
shape and suction slot in the airfoil might "pull the shock through" 
in a local region near the airfoil and thus effect a transition 
throv.gli the potential limit solution to operation on an upper branch 
of the basic curves, (fig. 3). 


30 


MCA TIT No. 1328 


APPENDIX B 


TBS OOI'C’EESSIBILITy COEKECTION RULE OF GREENE 


Hic comprc3sxl)ility correction rulo of Greene (reference 2) 
can be derived as follovs: For a fixed airfoil ordinate Y, the 

right-hand side of equation (7) may be witten for both conpressible 
and incompressible flows to yield the equation 


in which^ for a given airfoil^ both the compressible and incompress- 
ible curvature functions C and Cj_ satisfy the same boundary 
conditions. 

The Greene rule corresponds to the equation, equivalent to 
equation (18) of reference 2: 


For 7 = 1.4, the integral in equation (43) can be evaluated in 
closed form. The result is: 



(41) 



Jl 


(42) 


or 



(43) 


r-v 



pdv = 


i 1 


Jl 


Jl 



(44) 


vhere 


I'lo 


3 

(1+0.2 Mq2) 


G (Mo) = 


(45) 
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F (cp) = (cos 49 + 9 cos 29 + 23) + 59 


(46) 


The functions 
rospoctivel;^. 


sin 9o = 

Vi 

sin 9 = V 
G (Mq) and F ( 9 ) 


VO. 2 Mq 

+ 0.2 

sin 9 q 
are plotted 


(47) 


in figvires 22 and 23, 


Transition from equation (41) to (42) requires a certain corre- 
spondence between the compressible curvature function C and the 
incompressible curvature function C^. In order to see this corre- 
spondence, consider the integrands of equation (41 ) plotted against 
V and Vi, respectively (fig. 24(a)). A one -to -one ccrrcspondcnco 
between v and Vj_ can evidently always be established such that 
the elemental areas making up the integrals in equation (41) are 
equal, as indicated by the cross-hatched elements in figure 24(a); 
thus. 


pv-1 

Cv 


dv 


vj-1 

CiVi 


dVi 


(49) 


Aside from cases where a largo local supersonic region exists next 
to the airfoil, in which case it is possible for the curvature to 
have a maximum away from the airfoil, as indicated by the dotted 
line in figure 24(b), the compressible curvature function C and 
tho incompressible curvature fxinction both start from the same 

value Ca at tho airfoil and decrease monotonically (for chordwise 
stations near that of maximim velocity) to zero at y = 00 or 

V =■ Vi = 1. Eenco, a one-to-one correspondence between v and v^ 
can bo established in figure 24(b) such that C = at corresponding 

V and Vj^. If this one-to-one correspondence is the same as that 

by which equation (49) was obtained from (41 ), the curvature functions 
cancel out of equation (49), which can then bo integrated to yield 
equation (42). 


In general, the two coirrospondoncos just discussed are not the 
same. In this case the correspondence in v and v^ loading to 
equation (49) can bo regarded as retained and a compressible curva- 
ture function assumed that is obtained by this correspondence from 
the incompressible curvature function (indicated by tho dot-dash 
lino in figure 24(b)). This procedure again yields equation (42) 
from ( 49 ) by cancelation of C and Ci and integration. The com- 
pressible curvature function thus set up satisfies the boundary 
conditions and constitutes the essential approximation of Greene's 
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ru.le (in addition to the negD.ect of the curvatv.re of the potential 
lines Involved in equation (7)), 

The potential limit curve corresponding to Greene's rule is 
that for which 

dv 


or, by equation (49) 

pv = 1 


ESFERMCSS 

1. Stodola, A.; Steam and Gas Turbines. Vol. II. McGraw-Hill Book 

Co., Inc,, 1927, p. 992, 

2. Greene, Leonard Michael: The Attenuation Method for Compressible 

Flow Systems, Jour. Aero. Sci., vol. 12, no. 3, July 1345, 
pp. 329-338. 

3. Kaplan, Carl: The Flow of a Compressible Fluid Fast a Curved 

Surface. MCA ARE Wo. 3KD2, 1943. 

4. Kaplan, Carl: The Flow of a Compressible Fluid Past a Circular 

Arc Profile. NACA APJR Wo. L4Gis, 1944. 

5. Tsien, Hsue-Shcn: The "Limiting Line" in Mixed Subsonic and 

Supersonic Flow of Compressible Fluids. WACA TW Wo. 961, 1944. 

6. Glauert, H. : The Elements of Aerofoil and Airscrew Theory. 

Univ. Press (Cambridge), 1926. 

7. Oswatitsch, KI., and Wieghardt, K. : Zur Abschatzung dor kritischen 

Machschen Zahl. Zentralo fur wissenschaftlichos Berichtswesen 
der Luftfahrtforschung des Generalluftzeugmoisters (ZWB), Tech, 
Berichte, Bd. 10, Wr. 5, May 1943, pp. 157-159. 

8. (a) Taylor, G, I., and Mac coll, J. W. : Two-Dimensional Plow 

past a Curved Surface. Vol. II of Aerodynamic Theory, 
div. H, ch. IV, sec. 6, W. F. Durand, ed., Julius Springer 
(Berlin), 1943 (C.I.T, reprint), pp. 245-249. 

(b) Taylor, G. I., and Maccoll, J. ¥.: The Motion past Thin Air- 

foils. Op. cit., sec. 1, pp. 234-236. 


(50) 

(51) 


MCA TN Wo. 1328 


33 


9. von Kaniian, Th.: Compressibility Effects in Aerodynamics. Jour. 

Aero. Sci., vol. 8, no. 9, July 1941, pp. 337-556. 

10. Garrick, I. E., and Kaplan, Carl: On the Flow of a Compressible 

Fluid by the Hodograph Method. I - Unification and Extension 
of Present-Day Results. NACA ACE No. L4C24, 1944. 

11. Eingleb, Friednich: Exakte Losungen der Different ialgleichungen 

einer adiabatischen Gasstromung. Z.f.a.M.M., Bd 20, Heft 4, 
Aug. 1940, pp. 185-198. 

12. Meyer, Th,: Ueber zweidimensionale Bevegungsvorgange in einem 

Gas, das mit Ueberschallgeschwindigkeit stromt. Forechungs- 
arbeiten auf dom Gebiete des Ingenieur-wesens, Heft 62, VDI- 
Verlag G-.m.b.H. (Berlin), 1908, pp. 31-67. 

13. Taylor, G. I.: The Flow of Air at High Speeds past Ciurved 

Surfaces. R.&M. No. 1381, British A. E.C., 1930. 

14. Stack, Jolm: Comi^ressible Flows in Aeronautics. Jour. Aero. 

Sci., vol. 12, no. 2, April 1945, pp. 127-143; discussion, 
pp. 143-148. 


NACA TN no . 1328 


34 


TABLE I 

BASIC CALCULATIONS BY CURVATUBS METHOD FOR POSITIVE VELOCITY INCREMENTS 


Loge V 

0.04 

0.08 

0,12 

0.16 

0.20 

0.24 

0.28 

0.32 

0.36 

0.40 

0.44 

V - 1 


0.0408 

0.0833 

0.1275 

0.1735 

0,2214 

0.2712 

0.3231 

0.3771 

0.4333 

0. 4918 

0.5527 

Mo 

T1 













Curvature function. 

c = c^ 

'^y 

0 

1.0 

0.0404 

0.0816 

0.1237 

0.1666 

0.2104 

0.2552 

0.3008 

0.3474 

0,3960 

0.4435 

0.4931 


1.2 

.0452 

.0915 

.1380 

.1071 

.2366 

.2872 

.3390 

.3920 

.4462 

.5017 

.5505 


1.4 

.0523 

.1059 

.1609 

.2173 

.2752 

.33 45 

.3954 

. 4579 

. 5220 

.5878 

.6554 


1.6 

.0642 

.1302 

.1982 

.2681 

.3401 

.4143 

.4906 

.5692 

.6601 

.7335 

.8194 


1.8 

.0910 

.1850 

.2823 

.3829 

. 4869 

.5945 

.7067 

.8209 

.9400 

1 . 0633 

1.1908 


1.9 

.1287 

.2624 

.4010 

.5447 

.6938 

,8405 

1.0088 

1.1753 

1.3479 

1.5271 

1.7130 

.4 

1.0 

.0369 

.0744 

.1124 

.1509 

.1899 

.2294 

.2694 

.3099 

.3609 

.3923 

.4341 


1.2 

.0414 

.0834 

.1262 

.1697 

.2139 

.2588 

.3045 

.3509 

,3980 

.4450 

.4943 


1.4 

.0479 

.0966 

.1465 

.1974 

.2493 

.3023 

.3564 

.4116 

.4678 

.6252 

.5a37 


1.6 

.0687 

.1109 

.1007 

.2441 

.3091 

.3757 

,44 41 

.5142 

.5061 

.6598 

.7354 


1.8 

.0833 

.1692 

.2580 

.3496 

.4440 

.5416 

.6423 

,7463 

.8536 

.9644 

1.0788 


1.9 

.1179 

.2402 

.3669 

.4982 

.6343 

.7752 

.9212 

1.0725 

1.2294 

1.3919 

1,5604 

.5 

1.0 

.0348 

.0699 

.1055 

.1413 

.1774 

.2137 

.2503 

.2871 

.3240 

.3610 

.3981 


1.2 

.0390 

,0786 

.1106 

.1691 

.2001 

.2 416 

.2835 

.3259 

.3686 

.4117 

.4551 


1.4 

.0452 

.0911 

.1378 

.1963 

.2337 

.2828 

.3327 

.3834 

.43 49 

.4071 

.5401 


1.6 

.0554 

.1121 

.1702 

.2295 

,2902 

.3523 

.4159 

.4808 

.6472 

.6151 

.6844 


1.8 

,0737 

.1697 

.2433 

.3294 

.4181 

.6096 

.6039 

.7011 

.8012 

.9045 

1.0109 


1.9 

.1114 

,2260 

.3463 

,4701 

.5982 

.7300 

.8681 

1.010 

1.1576 

1.3101 

1 . 4680 

.6 

1.0 

.0321 

.0643 

,0965 

.1287 

,1609 

.1931 

.2251 

.2569 

.2083 

.3195 

,3501 


1.2 

.0369 

.0721 

,1085 

.1452 

,1820 

.2189 

.2558 

.2928 

.3290 

.3666 

.4032 


1.4 

.0415 

.0837 

.1263 

.1695 

.2131 

.2571 

.3015 

.3463 

.3914 

.4368 

.4825 


1.6 

.0511 

.1032 

,1563 

.2104 

,2655 

.5216 

,3788 

.4369 

,4960 

.5562 

.6173 


1.8 

.0726 

.1472 

.2240 

.3030 

.3042 

. 4676 

,6535 

.6418 

.7326 

.0259 

.9219 


1.9 

.1028 

.2093 

.3194 

.4332 

.5510 

.6727 

.7986 

.9289 

1,0636 

1.2029 

1.3471 

.7 

1.0 

,0286 

.0568 

,0846 

.1122 

.1392 

.1657 

.1915 

.2165 

.2405 

.2634 

.2860 


1.2 

.0319 

.0637 

,0955 

.1269 

.1501 

.1888 

.2192 

.2489 

,2779 

.3061 

.3333 


1.4 

.0370 

.07 41 

.1114 

.1487 

.1860 

.2232 

.2603 

,2972 

.3337 

.3699 

.4055 


1.6 

,0455 

.0915 

.1382 

.1856 

.2332 

.2814 

.3301 

.3792 

.4266 

.4784 

.5284 


1.8 

.0647 

.1310 

.1989 

.2686 

.3398 

.4128 

.4076 

.6642 

.6426 

.7229 

.0051 



.0917 

.1065 

.2843 

.3853 

.4895 

.5970 

.7080 

.8226 

.9409 

1.0631 

1.1892 

.8 

1.0 

.0237 

.0466 

,0686 

,0896 

.1093 

.1276 

.1442 

.1588 

.1711 

.1007 

.1870 


1.2 

.0267 

.0525 

,0888 

.1020 

,1253 

.1473 

.1679 

.1868 

.2038 

.2106 

.2307 


1.4 

.0308 

.0612 

,0912 

.1205 

.1490 

.1767 

.2033 

.2287 

.2527 

.2751 

.2955 


1.6 

,0380 

.0760 

.1140 

.1517 

.1093 

.2266 

.2634 

.2998 

.3355 

.3706 

.4041 


1.8 

.0541 

.1092 

.1657 

.2223 

.2002 

.3390 

.3987 

.4594 

.5209 

.5832 

.6464 


1.9 

.0769 

.1561 

.2375 

.3212 

.4072 

.4957 

.5066 

.6801 

.7763 

.0752 

,9769 

.85 

1.0 

.0206 

.0400 

.0500 

• 0743 

.0688 

.1009 

.1102 

.1162 

.1177 

.1135 

.1008 


1.2 

.0231 

.0452 

.0660 

.0854 

.1031 

.1187 

.1320 

.1423 

.1492 

.1510 

.1488 


1.4 

.0268 

.0529 

,0779 

.1018 

.1243 

.1452 

,1642 

.1810 

.1952 

.2064 

.2139 


1.6 

.0331 

.0659 

.0900 

.1296 

.1603 

.1901 

,2188 

,2 461 

.2719 

.2959 

.3179 


1.8 

.0474 

.0952 

.1436 

.1923 

.2414 

.2907 

.3404 

.3903 

.4403 

.4904 

.5405 


1.9 

.0674 

.1366 

.2073 

.2798 

.3541 

.4300 

.5070 

.5876 

.6691 

.7527 

.8302 

.9 

1.0 

.0167 

.0315 

.0441 

.0530 

.0600 

.0615 

.0557 

.0348 

.0476 

.0080 

.1251 


1.2 

.0187 

.0350 

.0507 

.0632 

.0725 

.0779 

.0779 

.0698 

.0453 

.0534 

.1016 


1.4 

.0219 

.0422 

.0600 

.0773 

.0918 

.1027 

.1097 

.1116 

.1067 

.0911 

.0514 


1.6 

.0270 

.0531 

.0779 

.1010 

.1224 

.1416 

.1604 

.1721 

.1821 

.1070 

.1879 


1.8 

.0390 

.0777 

.1162 

.1543 

.1919 

.2288 

.2651 

,3004 

.33 47 

.3670 

.3995 


1.9 

^0566 

.1121 

^1696 

^2279 

^2871 

^3472 

>4082 

.4700 

.5327 

.5962 

.6606 




Curvature function. C 













V/ 





0 

1.0 

.0402 

.0808 

.1218 

.1633 

.2052 

.2476 

.2903 

.3335 

.3772 

.4214 

“ . 46^6 

.5 

1.0 

.0347 

.0693 

.1039 

.1304 

.1729 

.2072 

.2413 

.2752 

• 3088 

.3421 

.3750 

.6 

1.0 

.0319 

.0636 

.0950 

.1261 

.1560 

.1870 

.2168 

.2460 

.2744 

.3021 

,3290 

.7 

1.0 

.0283 

.0561 

.0833 

.1098 

.1355 

.1603 

.1842 

.2068 

.2282 

.2401 

.2664 

.8 

1.0 

,0235 

.0460 

.0675 

.0876 

.1062 

.1231 

.1381 

.1608 

.1609 

.1679 

.1712 
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TABLE I - CONTINUED 

BASIC CALCULATIONS BY CURVATURE METHOD FOR POSITIVE VELOCITY INCREMENTS - CONTINUED 


Loge V 


0.40 

D.52 

3.56 

0.60 ( 

D.64 ( 

0.60 

0.72 

0.76 1 

0.80 

0.04 < 

3.88 

V - 1 


0.6161 

1.6020 

3.7507 

0.8221 ( 

D.0965 < 

0.9739 

1.0544 

1.1383 

1.2255 

1.3164 

1.4109 

Mo 






Curvature function, C 

- Ca(Y^) 





0 

1.0 

0.6430 

D.5955 

3.6484 

0.7023 ( 

0.7575 

0.0130 

0.0713 

0.9301 

0.9902 

1.0516 

1.1143 


1.2 

.6166 

.6761 

.7370 

.7994 

.0632 

.9206 

.9955 

1.064 

1.1343 

1.2062 

1.2798 


1.4 

.7247 

.7950 

.8609 

,9439 

1.0209 

1.0100 

1.1812 

1.2646 

1.3502 

1.4383 

1.5286 


1.6 

.9078 

.9990 

1.0929 

1.1896 

1.2893 

1.3920 

1.4978 

1.6070 

1.7194 

1.0354 

1.9549 


1.8 

1.3227 

1.4593 

1.6006 

1.7469 

1.8993 

2.0550 

2.2172 

2.3852 

2.5509 

2.7390 

2.9252 


1.9 

1.9058 

2,1059 

2.3136 

2,5291 

2.7527 

2,9040 

3.2256 

3.4758 

3.7350 

4.0045 

4.2839 

.4 

1.0 

.4764 

.5190 

. 5620 

.6053 

.6489 

.6927 

.7366 

.7808 

.8249 

.8691 

.9132 


1.2 

.5435 

.5933 

.6430 

.6949 

.7467 

.7989 

.0518 

,9052 

.9590 

1.0132 

1.0670 


1.4 

.6433 

.7041 

.7659 

.8269 

.8931 

.9583 

1.0247 

1.0922 

1.1607 

1.2305 

1.3012 


1.6 

.8129 

.092 4 

.9730 

1.0673 

1.1428 

1,2305 

1.3202 

1.4123 

1.5064 

1 . 6030 

1.7017 


1.8 

1.1968 

1.3107 

1,4446 

1,5746 

1,7087 

1.8472 

1.9900 

2,1378 

2.2901 

2.4476 

2.6099 


1.9 

1.7350 

1.9169 

2.1035 

2.2979 

2.4994 

2.7002 

2.9247 

3.1493 

3.3817 

3.6230 

3.0730 

.5 

1.0 

.4352 

.4722 

.5092 

.5459 

.5823 

.6105 

.6554 

.6070 

.7222 

.7559 

.7887 


1.2 

.4989 

.5420 

.5069 

.6312 

.6755 

.7199 

.7641 

.8003 

,8522 

.0959 

.9391 


1.4 

.5930 

.6482 

.7033 

.7690 

.8153 

.8722 

.9297 

.9076 

1.0460 

1.1048 

1.1640 


1.6 

.7553 

.8276 

.9016 

.9770 

1.0540 

1.1326 

1.2127 

1.2945 

1.3770 

1.4629 

1.5495 


1.8 

1.1206 

1.2337 

1.3602 

1,4703 

1,5941 

1.7216 

1.8531 

1.9007 

2.1283 

2.2724 

2.4207 


1.9 

1.6316 

1.8010 

1.9764 

2.1501 

2.3463 

2.5412 

2.7430 

2.9523 

3.1687 

3.3933 

3.6257 

.6 

1.0 

.3002 

.4097 

.4303 

.4661 

,4928 

.5183 

.5425 

.5652 

,5861 

.6052 

.6221 


1.2 

.4396 

.4756 

.5111 

.5461 

.5005 

.6140 

.6466 

,6783 

.7008 

.7380 

.7657 


1.4 

.5203 

.57 42 

.6203 

.6662 

.7121 

,7579 

.0034 

.0406 

.8934 

.9378 

.9017 


1.6 

.6793 

.7423 

.8063 

.0711 

.9369 

1.0035 

1.0710 

1.1395 

1,2086 

1.2787 

1.3496 


1.8 

1.0206 

1 . 1220 

1.2263 

1.3335 

1.4437 

1.5571 

1.6735 

1.7934 

1.9165 

2.0433 

2.1735 


1-9 

1.4962 

1.6504 

1.8100 

1.9751 

2.1459 

2.3226 

2.5054 

2.6946 

2.0904 

3.0932 

3 . 3030 

.7 

1.0 

.3237 

.3404 

.3540 

.3668 

,3759 

.3819 

.3042 

.3822 

.3754] 

.J629 

• 343S 


1.2 

,3594 

,3843 

.4076 

,4293 

.4491 

.4668 

.4021 

.4949 

.5047 

.5113 

• 51 45 


1.4 

.4406 

.4749 

.5084 

.5408 

.5722 

.6023 

.6310 

.6581 

.6836 

.7074 

.7292 


1.6 

.5787 

.6291 

.6797 

.7302 

,7800 

.8314 

.8017 

.9321 

.9821 

1.0321 

1.0019 


1.8 

,8893 

.9754 

1.0636 

1.1539 

1.2462 

1.3400 

1.4376 

1.5369 

1,6384 

1.7427 

1.8 495 


1.9 

1.3194 

1.4530 

1.5926 

1.7360 

1.8841 

2.0371 

2.1951 

2.3506 

2.5272 

2.7019 

2.8824 

.8 

1.0 

.1894 

.1860 

.178^ 

.1603 

.1205 

.0600 

.1160 

.1874 

.2484 

.3060 

.3619 


1.2 

,2398 

.2452 

.2462 

.2419 

.2309 

.2108 

.1771 

.1169 

.0957 

.1956 

.2676 


1.4 

.3139 

.3297 

.3427 

.3524 

,3504 

.3602 

.3569 

.3476 

.3313 

.3059 

.2684 


1.6 

.4376 

.4696 

.5000 

.5291 

.5565 

.6822 

.6059 

.6276 

.6472 

,6581 

.6797 


1.8 

.7105 

.7754 

.8411 

.9077 

.9752 

1.0436 

1.1129 

1.1834 

1.2549 

1.3278 

1 .4020 


1.9 

1.0815 

1,1891 

1,2998 

1,4138 

1,5311 

1.6519 

1,7764 

1.9047 

2.0371 

2.1739 

2.3151 

.85 

1.0 

.0725 

.0475 

.1135 

.1645“ 

.2129 

.260^ 

.3095 

.3509 

,4091 

.4601 

• 51 17^ 


1.2 

.1382 

.1158 

.0675 

,0895 

.1576 

,2154 

.2707 

.3255 

.3004 

.4357 

,4913 


1.4 

.2169 

.2145 

.2050 

.1858 

.1513 

.0810 

.1220 

.2063 

,2763 

.3417 

.4050 


1.6 

.3376 

.3547 

.3600 

.3795 

,3062 

.3086 

.3858 

.3772 

.3617 

,3377 

.3030 


1.8 

,5907 

.6408 

.6900 

.7406 

.7904 

.8401 

,6897 

,9393 

.9889 

1.0390 

1.0096 


1.9 

.9258 

1.0156 

1.3076 

1.2019 

1.2907 

1.3900 

1 . 5000 

1.6050 

1.7129 

1.8244 

1 . 93 93 

,g~~ 

■ ilo" 

.1628 

■ .20'4^ 

■ .24T^ 

.2670 

■ ,3312 

.3772 

, 42 46 

. 473 6 

,5240 

.5749 

. 6264 


1.2 

.1445 

.1098 

.2335 

.2790 

,3262 

. 3753 

.4261 

.4785 

,5322 

.5869 

, 6425 


1.4 

,0785 

.1402 

.1917 

.2426 

,2944 

.3476 

.4023 

.4586 

.5165 

.5751 

.6349 


1.6 

.1807 

.1605 

.1251 

.0281 

.1397 

.2144 

.2815 

.3463 

,4105 

.47 40 

. 5392 


1.0 

,4297 

.4571 

.4035 

.5070 

.5299 

.5496 

.5668 

,5816 

,5936 

.6034 

.6110 


1.9 

.7259 

.7914 

.8585 

.9266 

.9956 

1.0659 

1.1375 

1.2107 

1,2855 

1,3625 

1.4417 










NAT lONAL ADV ISORY 




Curvature function, u 

-^^Tc 


COMM 

ITTEE F 

OR AERC 

)NAUT ICS 

Q 

1.0 

.5112 

.5568 

. 6030 

.6497 

.6969 

.7447 

.7930 

.8420 

.8914 

.9415 

.9922 

i; 

1.0 

, 407 4 

.4393 

.4706 

,5012 

.5311 

.5600 

1 .5879 

.6147 

.6402 

. 6644 

.6870 

• 

A 

1.0 

.3540 

.3796 

. .4031 

.4253 

.4459 

.4648 

1 ,4018 

i .4967 

. 5092 

.5192 

.5263 

• V 

7 

1.0 

.2828 

.2972 

.3091 

.3185 

.3247 

.3275 

.3263 

,3204 

.3089 

' .2905 

.2631 

• » 
.8 

1.0 

.1699 

.1629 

1 .1480 

► .1207 

.0642 

.0989 

I .1659 

1 .2223 

.2751 

.5263 

.3766 
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TABLE I - CONCLUDED 


BASIC CALCUT,ATTONS BY CURVATUPE METHOD FOR POSITIVE VELOCITY INCREMENTS - CONCLUDED 


Loge V 
V - 1 

0.92 

1.5093 

0.96 

1.6117 

1.00 

1.7183 

1.04 

1.8292 

1.08 

1.9447 

1.12 

2.0649 

1.16 

2.1099 

1.20 

2.3201 

1.24 

2.4556 

1.28 

2.5966 

1.32 

2.7434 

Mq 

T\ 


{Yc; 





Curvature function, C = ^"0^ 

0 

1.0 

1.2 

1.4 

1.6 

1.8 

1,9 

1.1784 

1.3653 

1.6215 

2.0781 

3.1182 

4.5741 

1.2439 

1.4326 

1.7168 

2.2051 

3.3180 

4.8752 

1.3108 

1.5118 

1.0148 

2.3361 

3.5249 

5.1079 









.4 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

.9572 

1.1228 

1.3730 

1.8029 

2.7777 

4.1321 

1.0010 
1.1780 
1. 4458 
1.9064 
2 . 9508 
4.4008 

1.0444 

1.2334 

1.5197 

2.0123 

3.1296 

4.6794 

1.0875 

1.2889 

1.5945 

2.1208 

3.3141 

4,9680 

T 7 T 30 T 

1,3446 

1.6703 

2.2317 

3.5048 

5.2675 

TTTW 

1.4002 

1.7471 

2.3452 

3.7016 

5.5779 

1.2133 

1.4558 

1.8248 

2.4613 

3.9047 

5.0996 

■1.553 T 

1.5113 

1.9033 

2.5802 

4.1146 

6.2334 

1.2936 

1.5665 

1.9828 

2.7018 

4,3313 

6.5794 

1.3321 

1.6216 

2.0632 

2.0261 

4.5551 

6.9382 

1.369 C 

1.6763 

2.1444 

2.9535 

4.7864 

7,3104 

.5 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

.8205 

.9818 

1.2235 

1.6378 

2,5737 

3.8665 

.8511 

1.0241 

1.2832 

1.7278 

2.7314 

4.1160 

.8805 

1.0656 

1.3431 

1.8194 

2.8940 

4.3746 

,9086 

1.1064 

1.4033 

1.9128 

3.0616 

4.6424 

.9351 

1.1464 

1,4636 

2.0000 

3.2346 

4.9201 

.9600 

1.1854 

1.5241 

2.1051 

3.4128 

5.2079 

.9031 

1,2236 

1.5847 

2.2039 

3,5967 

5.506 C 

1,0045 

1.2600 

1.6455 

2.3048 

3.7866 

5.8152 

1.0240 

1.2970 

1.7065 

2,4078 

3.9826 

6.1360 

1,0416 

1.3323 

1.7677 

2.5131 

4.1649 

6.4681 

1.0576 
1.3667 
1.8294 
2 . 6207 
4.3940 
6.0129 

.6 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

.6368 

.7920 

1.0249 

1.4213 

2.3076 

3.5202 

.6491 
.8165 
1.0675 
1 . 4939 
2. 4456 
3.7451 

,6586 

.8394 

1.1095 

1.5675 

2.5876 

3.9781 

,6653 

.8603 

1.1508 

1.6420 

2.7339 

4.2192 

,6690 

.8795 

1.1914 

1.7176 

2.0847 

4.4693 

.6696 

.8960 

1.2315 

1.7944 

3.0401 

4.7284 

.6671 

.9123 

1.2712 

1.8725 

3.2004 

4.9968 

.6614 

.9261 

1.3105 

1.9523 

3,3660 

5.2753 

,6527 

.9386 

1.3499 

2.0337 

3.5371 

5.5642 

.6414 

.9501 

1.3895 

2.1172 

3.7139 

6.8638 

.6277 

.9610 

1.4290 

2.2031 

3.8971 

6.1750 

.7 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

.3435 

.5137 

.7491 

1.1316 

1.9592 

3.0692 

.3155 

.5088 

.7670 

1.1813 

2.0719 

3.2626 

.2759 

.4994 

.7831 

1.2310 

2.1078 

3.4628 

.2185 
.4051 
.79 72 
1.2011 
2.3071 
3.6701 

.1206 
. 4657 
.8098 
1,3316 
2.4302 
3.8851 

.1521 

.4407 

.0211 

1.3829 

2.5573 

4.1081 

.2544 

.4098 

.8315 

1.4353 

2.6086 

4.3392 

.3311 

,3725 

.8415 

1.4892 

2.8246 

4.5794 




.8 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

.4170 

.3341 

.2120 

.6929 

1.4781 

2.4611 

.4714 

.3970 

.1107 

.7043 

1.5561 

2.6123 

.5851 
.457 4 
.1609 
.7142 
1.6363 
2.7690 

.5770 

.5159 

.2629 

,7229 

1.7190 

2.9314 

.6293 
.5725 
.3419 
,7311 
1.8048 
3 . 1002 







.85 

1.0 

1.2 

1.4 

1.6 

1.8 

If® 

.5635 

.5472 

.4672 

.2530 

1.1410 

2.0582 

.6163 

.6031 

.5284 

.1756 

1.1936 

2.1814 

.6666 

.6586 

.5888 

.0859 

1.2478 

2.3092 









.9 

Tfo ^ 

1.2 

1.4 

1.6 

1.8 

1.9 

.6780 

.6983 

.7560 

.6038 

.6168 

1.5236 

.7293 

.7542 

.8126 

.6684 

.6214 

1.6085 
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Curvature function, C = C^/^— ^£-3. ^ 
\loge V2 

0 

.5 

.6 

,7 

.8 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0435 
.7078 
.5303 
.2230 
. .4262 

1.0954 
.7269 
. 5308 
.1604 
.4751 

1.1480 

.7439 

.5275 

.0534 

,5229 

.7536 

.5201 

.1694 

.5604 

.7710 

.5081 

.2702 

.6143 

.7006 

.4911 

.3374 

.7879 
. 4686 
.3970 

.7921 
,43 99 
, 4512 

.7932 

.4044 

.7913 

.3607 

.7862 

.3067 

Doge 
/ - 1 

V 

1.36 

2.8962 

1.40 

3.0552 

1. 44 
3.2207 

1.48 

3.3929 

1.52 

3.5722 

1.56 

3.7588 

1.60 

3.9530 

1.64 

4.1552 

1.60 

4,3656 

1.72 

4.6845 


Vo 

D 


Curvature function, C = Cg^^ ~ ■ 

.4 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

1.4061 

1,7308 

2.2265 

3.0838 

5.0253 

7.6965 

1.4413 

1.7848 

2.3096 

3.2172 

5.2721 

8.0968 

1.4753 

1.8386 

2.3937 

3.3539 

5.5273 

8.5121 

1.5081 

1,8920 

2.4788 

3.4939 

5.7909 

8.9429 

1.5398 

1,9453 

2,5652 

3.6375 

6.0636 

9-3899 

1.5705 
1 . 9985 
2 . 6528 
3.7849 
6.3458 
9.8537 

1 . 6003 
2.0518 
2.7420 
3.93 G 3 
6.6373 
10.3354 

1.6295 

2,1055 

2.8329 

4.0921 

6.9397 

10*8353 

1.6585 

2,1597 

2.9250 

4.2524 

7.2522 

11*3540 

1.6076 
2.2143 
3. 0210 
4. 4175 
7.5761 
11*8923 


.5 

1.0 

1.2 

1.4 

1.6 

1.8 

1.9 

1.0719 

1.4005 

1.8916 

2,7309 

4.6102 

7.1707 

1.0846 
1. 4338 
1.9546 
2.8430 
4.8336 
7.5417 

1.0964 

1.4668 

2.0186 

2.9599 

5.0648 

7.9269 

1.1074 

1.5000 

2.0838 

3.0792 

5,3040 

0,3265 

1.1181 

1.5337 

2.1507 

3,2022 

5.5518 

0.7416 







Curvature function, C * Ca^°-^ 

.5 

1.0 

0.7780 

0,7669 

0.7531 

0,7371 

0,7195 
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table II 

BASIC CALCULATIONS BY CURVATUBE METHOD FOR NEGATT\TE VIILOCITY INCREMENTS 


Loge V 
V - 1 

- 0.04 

0.0392 

- 0.12 

0.1131 

- 0.20 

0 . 1B13 

- 0.28 
0.2442 ■ 

- 0.36 
-a 3023 

- 0.44 

0.3560 

- 0.52 

- 0.4055 

- 0.60 

- 0.4512 

Mo 

T) 



Curvature function 

. c = c . 


r\ 

0 

1.0 

0.0396 

0.1165 

0.1904 

0.2615 

0.3299 

0.3958 

0.4592 

0.5203 


1.2 

.0443 

.1299 

.2118 

.2902 

.3654 

.4374 

.5065 

.5729 


1.4 

.0510 

• 1494 

.2429 

.3320 

.4169 

.4979 

.5751 

.6488 


1.6 

.0624 

.1820 

.2950 

.4018 

.5029 

.5986 

.6892 

.7751 


1.8 

.0880 

.2555 

.4124 

.5592 

.6967 

.8256 

.9465 

1.0599 


1.9 

.1243 

.3597 

.5787 

.7824 

.9721 

1.1406 

1.3131 

1.4663 

4.4 

1.0 

.0366 

.1077 

.1767 

.2437 

.3086 

.3716 

.4325 

.4915 


1.2 

.0406 

.1199 

.1963 

.2700 

.3411 

.4097 

.4759 

.5396 


1.4 

.0469 

.1377 

.2248 

.3082 

.3883 

.4650 

.5385 

.6090 


1.6 

• 0573 

.1676 

.2724 

.3720 

.4667 

.5568 

.6425 

.7241 


1.8 

.0808 

.2340 

.3796 

.5156 

.6434 

.7635 

.8666 

.9829 


1.9 

.1140 

. 3302 

.5317 

.7195 

.3947 

1.0580 

1.2105 

1.3528 

.5 

1.0 

.0345 

.1022 

.1684 

.2330 

.2958 

.3570 

.4165 

.4743 


1.2 

.0385 

.1138 

.1869 

.2578 

.3265 

.3930 

.4574 

.5197 


1.4 

.0444 

.1307 

.2138 

.2939 

.3710 

.4451 

.5165 

.5852 


1.6 

.0541 

.1588 

.2587 

.3540 

.4446 

.5316 

.6144 

.6934 


1.8 

.0764 

.2223 

.3598 

.4892 

.6112 

.7261 

.8344 

.9366 


1.9 

.1073 

.3123 

.5033 

.6815 

.8480 

1.0034 

1.1487 

1.2845 

.6 

1.0 

.0319 

.0953 

.1577 

.2191 

.2793 

.3382 

.3958 

.4521 


1.2 

.0356 

.1059 

.1748 

.2420 

.3076 

.3716 

.4337 

.4942 


1.4 

.0411 

.1215 

.1996 

.2753 

.3486 

.4196 

.4882 

.5545 


1.6 

.0501 

.1475 

.2409 

.3306 

.4166 

.4990 

.5780 

.6537 


1.8 

.0706 

.2060 

.3340 

.4550 

.5693 

.6774 

.7797 

.8764 


1.9 

.0996 

.2890 

.4662 

.6319 

.7870 

.9322 

1.0681 

1.1954 

.7 

1.0 

.0286 

• 0862 

.1438 

.2012 

.2581 

.3143 

.3697 

.4240 


1.2 

.0319 

.0958 

.1591 

.2217 

.2834 

.3441 

.4035 

.4617 


1.4 

.0367 

.1096 

.1812 

.2513 

.3199 

.3868 

.4620 

.5153 


1.6 

.0448 

.1327 

.2179 

.3004 

.3801 

.4571 

.5314 

.6031 


1.8 

.0631 

.1848 

.3005 

.4106 

.5152 

.6147 

.7092 

.7991 


1.9 

.0890 

.2587 

.4180 

.5675 

.7080 

.8400 

.9639 

1.0804 

• 8 

1.0 

.0243 

.0746 

.1261 

.1734 

.2312 

.2840 

.3366 

.3887 


1.2 

.0270 

.0825 

.1389 

.1957 

.2526 

.3092 

.3653 

.4207 


1.4 

.0312 

.0941 

.1574 

.2205 

.2831 

.3450 

.4059 

.4658 


1.6 

.0380 

.1135 

.1880 

.2613 

.3331 

.4033 

.4718 

.5384 


1.8 

.0532 

.1568 

.2566 

.3525 

.4447 

.5331 

.6179 

.6991 


1.9 

.0749 

.2186 

.3544 

.4828 

.6042 

.7190 

.8675 

.9301 

• 85 

1.0 

.0216 

.0669 

.1147 

.1641 

.2145 

.2654 

.3164 

.3673 


1.2 

.0240 

.0740 

.1260 

.1793 

.2334 

.2877 

.3419 

.3958 


1.4 

.0275 

.0841 

.1422 

.2010 

.2601 

.3186 

.3771 

.4349 


1.6 

.0335 

.1010 

.1689 

.2359 

.3030 

.3691 

.4341 

.4978 


1.3 

.0469 

.1389 

.2286 

.3156 

.4000 

.4817 

.5605 

.6366 


1.9 

.0659 

.1928 

.3137 

.4287 

.5382 

.6422 

.7412 

.8352 

• 9 

1.0 

.0182 

.0580 

.1015 

.1475 

.1952 

.2441 

.2936 

.3431 


1.2 

.0203 

.0638 

.1109 

.1603 

.2112 

.2630 

.3152 

.3675 


1.4 

.0230 

.0722 

.1242 

.1782 

.2334 

.2891 

.3450 

.4007 


1.6 

.0279 

.0860 

.1460 

.2072 

.2688 

.3305 

.3918 

.4524 


1.8 

.0390 

.1170 

.1946 

.2712 

.3466 

.4205 

.4927 

.5630 


1.9 

.0548 

.1612 

.2639 

.3628 

.4580 

.5494 

.6371 

.7212 
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Loge V 
V - 1 

- 0.08 

- 0.0769 

- 0.16 

- 0.1479 

- 0.24 

- 0.2134 

- 0.32 

- 0.2739 

- 0.40 

- 0.3297 

- 0.48 
- 0 . 3812 

- 0.56 

- 0.4288 

- 0.64 

- 0.4727 

Mo 



Curvature function, C = Cg 

/loge v\Ti 

vi 

0 

.5 

.6 

.7 

.8 

1.0 

1.0 

1.0 

1.0 

1.0 

0.0793 

.0693 

.0643 

.0680 

.0496 

0.1569 

.1382 

.1291 

.1172 

.1020 

0.2331 

.2068 

.1940 

.1775 

.1563 

0.3078 

.2750 

.2590 

.2384 

.2122 

0.3812 
• 3425 
.3237 
.2998 
.2692 

0.4533 
.4094 
.3882 
.3612 
• 3269 

0.5241 

.4757 

.4523 

.4226 

.3852 

0.5936 

.5412 

.5160 

.4839 

.4437 
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TABLE II - COHCLirDBD 

BASIC CALCULATIONS BY CURVATUNE METHOD POP NEGATIVE VELOCITY INCREMENTS - CONCLUDED 


cu 

cr> 

if \ 


Loge 
V - 1 


- 0,68 

- 0,4934 

- 0.76 

■ 0,6323 

- 0.84 

• 0.5683 

- 0.92 

- 0.6015 

- 1.00 

- 0.6321 

- 1.08 

■ 0.6604 

- 1.16 

- 0.6865 

- 0.7106 

Mo 




Curvature function , 

C = Ca 

(^r 

0 

1.0 

0.5792 

0.6361 

0.6909 

0,7439 

0.7951 

0.8445 

0.8924 

0,9387 


1.2 

,63 64 

.6976 

.7564 

.8129 

.0673 

.9197 

.9702 

1.0189 


1.4 

.7192 

.7864 

.0500 

.9123 

.9713 

1.0277 

1.0019 

1.1339 


1.6 

.0566 

.9339 

1.0073 

1.0770 

1.1434 

1.2065 

1.2666 

1.3239 


1.0 

1.1664 

1.2664 

1.3603 

1.4487 

1.5319 

1.6102 

1.6840 

1.7537 


1^9 

1.6090 

1.7422 

1.8663 

1.9823 

2.0905 

2.1915 

2.2860 

2.3744 

• 4 

1.0 

.5487 

.6041 

.6577 

.7097 

.7601 

.8089 

.8563 

.9022 


1.2 

.6012 

.6606 

.7177 

.7730 

.6264 

.8779 

.9277 

.9758 


1.4 

.6767 

.7416 

.0039 

.0638 

.9212 

.9765 

1.0296 

1.0807 


1.6 

,8017 

.8757 

.9462 

1.0134 

1.0775 

1.1387 

1.1972 

1.2530 


1.8 

1.0830 

1.1774 

1.2662 

1.3501 

1 . 4292 

1.5039 

1.5745 

1.6412 


1.9 

1.4057 

1.6099 

1.7259 

1.8344 

1.9350 

2.0308 

2.1198 

2.2031 

.5 

1.0 

.5304 

.5860 

.6379 

.6894 

.7393 

.7877 

.83 48 

.8805 


1.2 

.5800 

.6303 

.6947 

.7492 

.8020 

.0530 

.9024 

.9502 


1.4 

.6512 

.7148 

.7759 

.0340 

.8914 

.9459 

.9985 

1.0491 


1.6 

.7608 

.8408 

.9095 

.9753 

1.0301 

1.0982 

1.1557 

1.2100 


1.8 

1.0329 

1.1239 

1.2097 

1.2909 

1.3676 

1.4402 

1.5009 

1.5739 


1.9 

1.4114 

1.6302 

1.6413 

1.7454 

1.6428 

1.9342 

2.0196 

2.1002 

.6 

1.0 

.5070 

.5606 

.6126 

.6633 

.7127 

.7607 

.8075 

.0530 


1.2 

.5629 

.6099 

.6652 

.7180 

.7708 

.0213 

.8702 

.9176 


1.4 

.6185 

.6803 

.7400 

.7976 

.8532 

.9069 

.9507 

1.0080 


1.6 

.7263 

.7968 

.8624 

.9263 

.9875 

1.0463 

1.1026 

1.1567 


1.8 

.9680 

1.0546 

1.1366 

1.2144 

1.2081 

1.3579 

1.4242 

1.4872 


1.9 

1.3147 

1.4265 

1.5314 

1.6298 

1.7221 

1.0008 

1.8903 

1.9669 

.7 

1.6 

. 4*^74 

.^^96 

.5807 

.6306 

.6793 

.7269 

.7733 

• 8185 


1.2 

.5185 

.5739 

.6279 

.6806 

.7316 

.7814 

• 8298 

.0768 


1.4 

.5769 

.6366 

.6945 

.7606 

.0050 

.0577 

.9007 

.9581 


1.6 

.6721 

.7386 

.8026 

.8642 

.9234 

.9806 

1.0355 

1.0684 


1.0 

.8845 

.9656 

1.0430 

1.1166 

1.1865 

1.2530 

1.3165 

1.3769 


1.9 

1.1899 

1.2929 

1.3898 

1.4810 

1.5668 

1.6477 

1.7240 

1.7959 

• 8 

1.0 

.4402 

.4910 

.5409 

.5900 



"■7511 

.7761 


1.2 

.4753 

.6288 

.6813 

.6327 

.6830 

.7320 

.7798 

.8265 


1.4 

.5243 

.5816 

.6375 

.6919 

.7449 

.7965 

.8466 

• 6954 


1.6 

.6032 

.6660 

.7269 

.7859 

.8430 

.8982 

.9516 

1.0033 


1.0 

.7770 

.8518 

.9229 

.9913 

1.0567 

1.1194 

1.1794 

1.2369 


1.9 

1.0271 

1.1189 

1.2067 

1.2078 

1.3656 

1.4393 

1.5091 

1.5753 

• 85 

1.0 

.4178 

.4679 

.6172 

.5658 

.6135 

.6603 

.7062 

.7512 


1.2 

.4491 

.5017 

• 5534 

• 6042 

.6539 

- .7026 

.7502 

.7967 


1.4 

.4918 

.5477 

.6025 

.6561 

.7004 

.7594 

.8091 

.8575 


1.6 

.5601 

.6209 

.6800 

.7376 

.7935 

.8477 

.9003 

.9514 


1.0 

.7100 

.7307 

.8486 

.9141 

.9769 

1.0374 

1.0955 

1.1514 


If® 

.9246 

1.0095 

1.0903 

1.1671 

1.2401 

1.3095 

1.3755 

1.4384 

• 9 

1.0 

.3926 

.4419 

.4906 

• 5388 

.5862 

.6328 

• ^786 

.7236 


1.2 

.4195 

.4711 

.5220 

.5722 

.6215 

.6699 

.7173 

.7637 


1.4 

.4569 

.6104 

• 5641 

.6168 

.6684 

.7189 

.7682 

.6164 


1.6 

.5121 

.6707 

• 6281 

.6842 

.7309 

.7922 

.8441 

.6946 


1.0 

.6314 

.6977 

.7621 

.82 43 

.8845 

.9427 

.9989 

1.0632 


JLlL. 

.8016 

.8790 

.9629 

1.0237 

1.0913 

1.1561 

1.2100 

1.2774 
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Log, V 
V - 1 

-0.72 

-0.5133 

-0.80 

-0.5607 

-0.88 

-0.5852 

-0.96 

-0.6171 





Mo 



1 log* v\^ 

Curvature function, C ■ 7/ 

0 

• 6 
.6 

• • . • . 

0000 0 
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.6193 
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Figure !•- Compressible potential flow past symmetrical airfoil# 
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Airfoil velocity Increment, V-1 
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Airfoil curvature parameter, jYCa 


Figure 2.- Positive velocity Increment on airfoil as a function of airfoil ordinate 
y and airfoil curvature for various values of free-stream Mach number Mq. 
Values computed from equation (10) and tabulated In table !• r\ - 1,0, 
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(a) Pree-stream Mach nxunber Mq • 0. 

Plgtire 3.- Positive velocity Increment on airfoil as a function of airfoil ordinate Y and 
airfoil curvature Ca for various values of t). Values computed from equation (10) for 
V>1 and tabulated In table !• 




Airfoil velocity Increment* V-1 


F i g. 3b 


NACA TN no. 1328 



NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 




|.^4„ I.. 


Airfoil curvature parameter* 


1 


|:1'. 

t 

1 

f 

i 

1 

f- 

± 

til 


(b) Pree-stream Mach number M© * 0»4* 

Figure 3,- Continued, Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y^and airfoil curvature for various values of Values computed from 
equation (10) for V>1 and tabulated in table 1, 
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Airfoil curvature parameter, /y^ 
(o) Free -stream Mach number M^ * 0«5« 


Figure 3.- Continued* Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature C« for various values of Values computed from 

equation (10) for V >1 and tabulated in table I. 
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Fig, 3d 
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(d) Pree-stream Mach number * 0.6. 

Figure 3.- Continued. Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature Ca for various values of r\. Values computed from 
equation (10) for V>1 and tabulated In table I, 
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Figure 3.- Continued, Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature for various values of t^. Values computed from 
equation (10) for V>1 and tabulated In table I, 
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(f) Pree-stream Mach number Mq • 0.8. 

Flffure 3.- Continued. Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature Ca for various values of Values computed from 
equation (10) for V >1 and tabulated In table I, 
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(g) Pree-stream Mach number Mq « 0.85, 

Figure 3.- Continued. Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature C« for varIou 5 values of t). Values computed from 
equation (10) for V >1 and tabulated In tabic I, 
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Fig. 3h 
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Figure 3,- Concluded. Positive velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature Ca lor various values of n. Values computed from 
equation (10) for V>1 and tabulated in table I, 
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(a) Pree-stroaa Maeb nunber Mg * 0. 


Figure 4.» VeMtlve reloelty IncreMnt on airfoil at a function of airfoil ordinate 
I and airfoil eurratur# for rarlous valuea of n. Values computed from equation 
(10) for VCl and tabulated In table II. *rom equowon 


Fig. 4b 
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(b) Pree-stream Mach number Mq • 0.4. 

Figure 4.- Continued. Negative velocity Increment on airfoil as a function of 
airfoil ordinate Y and airfoil curvature for various values of *n. Values 
computed from equation (10) for V<1 and tabulated In table II. 
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(c) Pree-stream Mach number * 0.6. 


Figure 4.- Continued. Negative velocity Increment on airfoil as a function of airfoil 
Ordinate Y and airfoil curvature Cg for varloui values of Values computed from 
equation (10) for V<1 and tabulated In table II, 


J 


NACA TN NO- 1328 


Fig. 4d 


NATIONAL ADVISORY 
:COMMITTEC FOR AERONAUTICS- 



o 

U 


•3 *4 *5 »6 

Airfoil curvature parameter, -|YCa 


(d) Pree-streara Mach number Mq ■ 0.6. 


Figure 4.- Continued. Negative velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature for various values of ti. Values computed from 
equation (10) for V<1 and tabulated In table II, 
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Airfoil curvature parameter, 

(e) Pree-stream Mach number Mq ■ 0.7. 

Figure 4.- Continued* Negative velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature for various values of t). Values computed from 
equation (10) for V<1 and tabulated In table II, 
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(f) Free -stream Mach number Mq * 0.8. 

Figure 4.- Continued. Negative velocity Increment on airfoil as a function of airfoil 
ordinate T and airfoil curvature Ca for various values of -n. Values computed from 
equation (10) for V<1 and tabulated In table II. 
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(g) Pree-stream Mach number * 0.85, 

Figure 4,- Continued, Negative velocity Increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature for various values of t). Values computed from 
equation (10) for V<1 and tabulated In table II, 
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(h) Pree-stream Mach number M© » 0.9. 

Figure 4.- Concluded. Negative velocity increment on airfoil as a function of airfoil 
ordinate Y and airfoil curvature ^or various values of Values computed from 
equation (10) for V<1 and tabulated in table II. 
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Figure 5.- Basic data for Kaplan symmetrical section of 
thickness ratio 0. IQ (reference 3). 
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Fig. 6 
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Figure 6*- Velocity distributions on Kaplan section of 
thickness ratio 0,10 computed by equation (14), 
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(a) Relation between |[Y^ and r\. 

Figure 7.- Potential limit curves; loci of points of 
infinite slope in figure 3. 


Airfoil velocity increment, V-1 


Fig. 7b 
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(b) Relation between V-l and tj. 

Figure 7.- Continued. Potentiail limit curves; loci of 
points of infinite slope in figure 3. 
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F’ree-stream Mach number, 

(c) Relation between and Mq. 

Figure 7.- Continued. Potential limit curves; loci of 
points of infinite slope in figure 3. 
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Fig. 7d 
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Figure 7.- Concluded. Potential limit curves; loci of 
points of infinite slope in figure 3. 
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Chordwise losation, x, semichord 

(a) Mq = 0.5. 

Figure 8.- Comparison of velocity distribution on 
Kaplan section of thickness ratio 0. 10. 
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Chordwise location, x, semichord 
(b) Mq = 0.75. 

Figure 8.- Continued. Comparison of velocity distribution 
on Kaplan section of thickness ratio 0. 10. 
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Fig. 8c 



Chordwise location, x, semichord 
(c) Mq = 0.83. 

Figure 8.- Continued. Comparison of velocity distribution 
on Kaplan section of thickness ratio 0. 10. 
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Fig, 8d 
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(d) Mq = 0.843. 

Figure 8.- Concluded. Comparison of velocity distribution 
on Kaplan section of thictcness ratio 0. 10. 


V71 

VO 


Section ordinate 


NACA TN No. 1528 
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Figure 9.- Basic data for symmetriceil biconvex 
section of thickness ratio 0.20. 
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Chordwise location, x, semichord 

Figure 10.- Velocity distributions biconvex section 
of thickness ratio 0.20. 
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(a) Mo = 0. 


Figure II.- Streamline-curvature variation in field of airfoil. Kaplan section 

of thickness ratio 0. 10. 
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(b) Mq = 0.5. 


Figure II.- Continued. Streamline-curvature variation in field of airfoil. Kaplan 

section of thickness ratio 0. 10. 
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(c) Mq = 0.843. f 

Figure II.- Concluded. Streamline-curvature variation in field of airfoil. Kaplan _ 

section of thickness ratio 0. 10. o 
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Figure I 2 .- Velocity variation i n field of a i rfoil. Kaplan section ojf thickness 

ratio 0. 10. 
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Nondimensional distance normal to airfoil, ~ ^ 

(b) Mq = 0.5. 

Figure 12.- Continued. Velocity variation in field of airfoil. Kaplan section of 

thickness ratio 0. 10. 
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(c) Mq =• 0.843. 

Figure 12.- Concluded. Velocity variation in field of airfoil. Kaplan section of 

thickness ratio 0.10. 
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B'igure 14.- Local supersonic regions for Kaplan 10- 
percent section. 
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Figure 15,- Dependence of maximum velocity increment on 
circular arc camber line on position of reference point 
along airfoil streamline. The quantity f is the ratio 
in incompressible flow of the velocity decrement at a 
point on the airfoil streamline to the maximum velocity 
Increment on the circular arc. 


g- 


16 


NACA TN NO, 


528 



P 


u 

0) 

4^ 


a; 

B 




c 

o o 

«H4> 

•p 10 
O g 


Q>«H 

«4 

o 


Figure 16. - Basic data circular arc-mean camber line 
at zero angle of attack canioer ratio, 0.04. 
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Figure 17.- Comparison of velocity distributions on 
circular arc mean camber line by different methods. 
Camber ratio, 0.04, 
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Pressure coefficient 


Fig, 18 a 
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Free-streaflj Mach number, Mq 

(a) Negative-pressure increments. 

Figure 18.- Compressibility-correction rules by curvature method. Potential limit curves from 

figure 7(d). 
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Fig. 


Free-stream Mach number, Mq 

(b) Positive-pressure increments. 

Figure 18.- Concluded. Compressibility-correction rules by curvature method. Potential limit 

curves from figure 7(d). 
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Figure 19.- Comparison of compressibility correction rules. 



Figure 20.- One-dimensional curves for converging-diverging channel analogous to curves of figure 2. 

A = I /pV. 
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Figure 21,- Converging-diverging channel. 
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Figure 23,- Plot of function 

F(<p) * (cos 4«p 9 cos 2? 23) 5»- 
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(a) fdv = fidvi. 


(b) Curvature functions. 


Figure 24«- Correspondence between compressible and Incompressible flow for 
compressibility correction rule of Greene, 

4 ^ 
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